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Abstract 

We review the status of searches for Supersymmetry at the Tevatron 
CoUider. After discussing the theoretical aspects relevant to the produc- 
tion and decay of supersymmetric particles at the Tevatron, we present 
the current results for Runs la and lb as of the summer of 1997. 
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1 Introduction 



The Tevatron is a 1-kilometer radius superconducting accelerator ring, located 
at the Fermi National Accelerator Laboratory (Fermilab), in Batavia, Illinois, 
U.S.A. The ring is used in two modes: as a source of high energy beams for 
fixed target experiments, and, in conjunction with the Antiproton Source, as 
a proton-antiproton collider, operating with a pp center-of-mass energy of 
•^i = 1.8 TeV. The Tevatron Collider has been the world's accelerator-based 
high energy frontier since it first began taking data in 1987, and has thus been 
a prime location to search for the final pieces of the Standard Modelll (SM) 
and new phenomena beyond. _ 

With the discovery of the top quark at the Tevatronp the SM particle 
spectrum is almost complete, with only the Higgs boson (and, arguably, the 
tau neutrino) lacking direct experimental confirmation. If the interactions of 
the leptons, quarks and gauge bosons of the SM remain perturbative up to 
very high energies (as appears to be the case from the measured running of the 
gauge couplings), then the mechanism responsible for electroweak symmetry 
breaking (EWSB) is expected to contain one or more fundamental scalar Higgs 
bosons that are light, i.e. with masses of the order of the symmetry-breaking 
scale. 1^ 

The Higgs mechanisnJj plays a crucial role in the SM. The neutral compo- 
nent of the Higgs boson acquires a vacuum expectation value to give mass to 
the W and Z gauge bosons as well as to the SM fermions. In addition, the cou- 
plings of the Higgs boson to the gauge bosons and fermions are such as to cancel 
the infinities in the electroweak radiative corrections and prevent unitarity vio- 
lation in the longitudinal scattering of the gauge bosons. However, if the SM is 
vahd up to an energy cutoff of Acutoff « A'/pianck (Mpianck = 1-22 x 10^^ ^1^^' 
the model has a fundamental problem, the so-called naturalness prohlenv^ 

The radiative corrections to the Higgs boson mass-squared calculated in 
the SM are quadratically divergent (proportional to A^^^^g). A physical Higgs 
mass of the order of the electroweak scale requires a cancellation of one part 
in 10"'^^ between these radiative corrections, which come from the interactions 
of the Higgs bosons with all other particles in the theory, and the bare Higgs 
mass at the Planck scale: 

m]j ~ (Acutoff ) - a Acutoff ■ ( 1 ) 

Either there is an extreme "fine tuning" necessary to have a cancellation of 
two independent effects (the naturalness problem) , or there must be some new 
principle at work. nnn 

Supersymmetry (SUSY)Ho'Q is a new symmetry which provides a well— 
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motivated extension of the SM with an elegant solution to the naturalness 
problem. Supersymmetric transformations relate fermionic and bosonic de- 
grees of freedom. Each left-handed and right-handed fermion of the SM is 
postulated to have its own bosonic superpartner with equal mass and coupling 
strengths. Similarly each SM boson would have its own fermionic superpart- 
ner, again with equal mass and couplings. Because bosons and fermions induce 
radiative effects of opposite signs, SUSY naturally provides an exact cancella- 
tion of the otherwise quadratically-divergent radiative corrections to the Higgs 
boson massp 

Given that no superparticles have been observed so far, it is assumed that 
SUSY is broken, and that in general the sparticles must be heavier than their 
partners. In order to break SUSY without spoiling the necessary cancellation 
of quadratic divergences, the splittings between the masses of the SM particles 
and their SUSY partners should not be much larger than a few TeV. If SUSY 
is a consistent description of Nature, then the lower range of sparticle masses 
can be within the reach of thciITevatronp motivating a wide range of searches 
in a large number of channelsEj The jnass of the lightest neutral Higgs boson 
is strongly constraipjtLwithin SUSY^j and could be within the reach of the 
upgraded TevatronOEj In addition, the decay of the heavy top quark, pair 
produced in strong interactions with a cross section of ~ 6 pb at the Tevatron, 
gives a unique mechanism for producing lighter supersymmetric particles which 
might not otherwise be produced at a large rate in proton-antiproton collisions. 

Two experimental collaborations, CDF and D0, have large, general-purpose 
detectors at Fermilab. The Tevatron had initial running periods in 1985 and 
1987 with low luminosity.f] In 1988 - 1989 (the '89 Run'), the Tevatron oper- 
ated at \/s = 1.8 TeV with an average instantaneous luminosity of 1.6 x 10"^° 
cm~^s~^, and the CDF detector collected approximately 4.4 pb~^ of data. 
In 1992 — 93 (Run la), CDF and D0 accumulated approximately 20 and 15 
pb"^, respectively, and in 1994 — 95 (Run lb) 90 and 108 pb~^, for a total 
Run I integrated luminosity of more than 100 pb"^ per detector^] The average 
instantaneous luminosity during Run I was approximately 1 x 10^^ cm~^s~^. 
A new run (Run II) utilizing the new Main Injector and Recycler rings and 
other major accelerator improvements is scheduled to begin around the year 
2000, reaching an average instantaneous luminosity of 1-X lO'^^ cm~^s^^ and 
an expected integrated luminosity of 1000 pb~^ per yearta In addition, the en- 

"The 1985 run produced the first detected luminosity, with 20 events recorded by CDF. 

''The difference in integrated luminosities in Run lb comes partly from the fact that 
D0 uses the 11.3 pb~^ from Run Ic while CDF ignores it, and the fact that the D0 experi- 
ment normalizes its luminosity (and hence all cross sections) to an inelastic cross section that 
is 2.4% smaller than that used by CDF. For the actual luminosities used in each analysis see 
Table |. 
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ergy of the machine wiU be increased to -/i = 2 TcV, substantiaUy increasing 
the cross section for producing heavy particles (the top quark pair production 
cross section, for example, increases by 40% from y/s = 1.8 TeV to 2 TeV). 
Run II, with an upgraded Collider and detectors, holds a great deal of promise 
for Higgs boson and sparticle searches. 

2 The MSSM 

In the past two decades, a detailed picture of the Minixtml-Supersymmetric 
extension of the Standard Model (MSSM), has emerged.B0O In the MSSM, 
the particle spectrum is doubled by SUSY. Moreover, to generate masses for 
up- and down-type fermions while presepiing SUSY and gauge invariance, the 
Higgs sector must contain two doubletsEa After EWSB, there is a quintet of 
physical Higgs boson states: two CP-even scalar {h,W), one CP-odd pseu- 
doscalar (A), and a pair of charged (H"^) Higgs bosonsEj All the Higgs bosons 
and other SM particles have superpartners with the same quantum numbers 
under the SM gauge groups SU{3)c x SU{2)l x U{1)y, but with different 
spinel The spin-1/2 partners of the gauge bosons (gauginos) are denoted as 
winos , zinos Z, photinos'[| 7, and gluinos g. The spin-1/2 partners of the 
Higgs bosons (Higgsinos) are Hi,H2 and H^. Because of EWSB, the Hig- 
gsinos and SU{2)l x U{1)y gauginos mix to give physical mass eigenstates 
consisting of two Dirac fermions of electric charge one, the charginos xt2' ^^'^ 
four neutral Majorana fermions, the neutralinos Xi-4- The spin-0 partners of 
the fermions (sfermions) ^ are squarks Q, | sleptons £ and sneutrinos v. Each 
charged lepton or quark has two scalar partners, one associated with each chi- 
rality. These are named left-handed squarks and sleptons, which belong to 
SU{2)l doublets, and right-handed squarks and sleptons, which are SU{2)l 
singlets. The neutrinos have only left-handed superpartners v, which belong 
to SU{2)l doublets. The gluino g and squarks Q carry color indices and are 
SU{3)c octets and triplets, respectively. 

The MSSM Lagrangian contains interactions between particles and spar- 
ticles, fixed by SUSY. There are also a number of soft SUSY-breaking mass 
parameters. "Soft" means that they break the mass degeneracy between SM 
particles and their SUSY partners without reintroducing quadratic divergences 

"The superpartners of the U{1)y and SU{2)i^ gauge bosons (before EWSB) are the Bino 
B, the unmixed neutral Wino W3, and the unmixed charged Winos Wi and W2- 
''Charge conjugate scalars are denoted by *, e.g. Q* . 

"^To allow easier reading, we use the non-standard symbol Q instead of q. This has no 
special significance. 
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while respecting the gauge invariance of the theory. The soft SUSY-breaking 
parameters are extra mass terms for gauginos and scalar fermions, and trilinear 
scalar couplings. The exact number of extra parameters depends on the exact 
mechanism of SUSY breaking. In the remainder of this section, the MSSM 
particle spectrum and properties will be described in general, and also with 
reference to specific SUSY-breaking scenarios and variations. 



2.1 Sparticle Spectrum 

The chargino and ncutralino masses and their mixing angles (that is, their 
gaugino and Higgsino composition) are determined by the SM gauge boson 
masses {Mw and Mz), tan/3,|]two soft SUSY-breaking parameters (the SU{2)l 
gaugino mass M2 and the U{1)y gaugino mass Mi), and the SUSY Higgsino 
mass parameter /i, all evaluated at the electroweak scale Af^jvy-f] Explicit so- 
lutions are found by considering the 2x2 chargino Mc and 4x4 neutralino 
Mn mass matrices: Q 

^ M2 V2Mws(3 

^ \ V2Mwcp M 

/Afi 0^M=^^° '^'^ 7=f -^zcpsw Mzs(3sw \ 
' \ Q M2 J' ^' \ -fi )■■ \ Mzcficw -MzsPcw J 

Mc is written in the W+ - H+ basis, Mn in the B-W^-Hi-H2 basis, with 
the notation s/3 = sin /?, c/3 = cos /3, sw — sin 9w and cw = cos 9w- In general, 
the mass eigenstates are admixtures of the interaction states, but, for large 
values of or Mi and A/2, the limit is reached where the mass eigenstates are 
mostly pure gaugino or Higgsino states (independent of tan/3). In particular, 
if 1^1 ^ Mz and Afi,Af2 — Mz, with Mi < M2, the lightest eigenstates are 
gaugino-like and the heaviest are Higgsino-like, leading to the spectrum: 

Af~± ~ A/2; M~± ~ \^^\ (3) 
M~o ^ Mi; M~o ^ M2; M~o ^ M~o ^ ImI- 

Ai A2 A3 A4 

'*One Higgs doublet, H2 , couples to it, c, and t, while the other. Hi , couples to d, s,b,e, /x, 
and T. The parameter tan/3 is the ratio of vacuum expectation values {H2) / {Hi) = V2/V1, 
and v'^ = v'^ + v^, where v is the order parameter of EWSB. 

''The electroweak scale Mew is roughly the scale of the sparticle masses themselves. 
Usually, in the literature, for simplicity, M^w — Mz- 

^Beware of different sign conventions for fj, in the literature. Both PYTHIA and ISA JET use 
the convention used here. 
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Similarly, if Afi, M2 ^ Mz and \n\ — Mz, the lightest eigenstates are Higgsino- 
like and the heaviest are gaugino-like: 

M~± ~ ImI; M~± ^ A/2 (4) 
M~o ~ M~o ^ M~o ~ Mi; M~o ~ M2. 

Ai A2 A3 A4 

Another interesting example where Xi is Higgsino-like, X2 is photino-like, and 
all other charginos and neutralinos are mixtures, occurs for Mi = A/2 ~ ~ 
Mz and tan/3 ~ 1: 



A/~±^ = 1 1 A/2 + /i T x/ (M2 - + 4A/2, 1 (5) 



A/~o = |/.|; A/~o=A/i; A/~o = i|A/2 + ^ T J {M^ - m)^ + 4A/||. 

Al A2 A3_4 Z » 

Since the SU {3)c symmetry of the SM is not broken, the gluinos have masses 
determined by the SU(3)c gaugino mass parameter A/3I] The neutralinos and 
the gluinos are Majorana particles, and do not distinguish between states and 
their charged conjugate. Depending on their Higgsino or gaugino composition, 
the X couplings to gauge bosons, and left- and right-handed sfermions will dif- 
fer substantially, and the production and decay processes will strongly depend 
on that composition (see the discussion below). 

The mass eigenstates of squarks and sleptons are, in principle, mixtures of 
their left- and right-handed components, given for the first generation by: 

™ii - "iQi + Du^ ml^ ~ mfj^ + ml + Du^ 

"lie - "^ii + -^Se> (6) 

where Jtiq^ , to|^^ , m^^ , m'jj_^ , and m|;^ are soft SUSY-breaking parameters 
and Df^ = A/| cos(2/3)(T3^ - Qfsiii^0w), = A/| cos(2/3)Q/ sin^ 6*^/ are 
D-terms Q associated with EWSB {T^^, is the weak isospin eigenvalue of the 
fermion, Qf the electric charge). A similar expression holds for the second 
(third) generation with the substitutions u c{t),d s{b),e — > /i(T),l —^ 

^The physical gluino mass is shifted from the value of the gluino mass parameter M3 
because of radiative corrections. As a result, there is an indirect dependence on the squark 
masses. 

''D— terms are terms in the scalar potential which arc quartic in the fields and are pro- 
portional to the gauge couplings squared. 



6 



2(3). In most high-energy models, the soft SUSY-breaking sfermion mass pa- 
rameters are taken to be equal at the high-energy scale, but, in principle, they 
can be different for each generation or even within a generation. However, the 
sfermion flavor dependence can have important effects on low-energy observ- 
ables, and it is often strongly constrained. The suppression of flavor changing 
neutral currents (FCNC's), such as Kl Tr°vD, requires that either (i) the 
squark soft-SUSY breaking mass matrix is diagonal and degenerate, oil (ii) 
the masses of the first- and second-generation sfermions are very large. EEI 

The magnitude of left-right sfermion mixing is always proportional to 
the mass of the corresponding fermion. The left-right mixing of squarks and 
sleptons of the first and second generation is thus negligible, and Ql.r, with 
Q = u,d,c,s, and iL,R,i'i, with £ = e,/Lt, are the real mass eigenstates with 
masses m~ and m^^ ^, m^j^,, respectively. For the third generation sfermions, 
the left-right mixing can be nontrivial. The mass matrix for the top squarks 
(stops) in the {iL,iji) basis is given by 

= f + ™? + A\ rntiAt ~ ^i/ tan/3) \ 

''^t mt (At - m/ tan (3) +mi + D^^ ^' > 

where At is a soft SUSY-breaking parameter| Unless there is a cancellation 
between At and /i/tan/3, left-right mixing occurs for the stop squarks because 
of the large top quark mass. The stop mass eigenstates are then given by 

ti = cos tL+sin^f in 

t2 = — sin0( Il+cosOi Ir, (8) 

where the masses and mixing angle 9^ are fixed by diagonalizing the squared- 
mass matrix Eq. (|^). Because of the large mixing, the lightest stop ii can be 
one of the lightest sparticles. For the sbottom, an analogous formula for the 
mass matrix holds with mjj^ — > ^Di, At Ajj^ D^^ ^ — s- _D^^ ^, rrit mf,, 
and tan/3 — s- l/tan/3. For the stau, substitute mg^ rriLs, "^u^ ~^ i^Ea, 
At — > _Dj^ ^ — !■ Dfj^ ^, mt — > rjir and tan/3 l/tan/3. The parameters 
At, Af,, and Ar can be independent soft SUSY-breaking parameters, or they 
might be related by some underlying principle. When mi, tan (3 or tan /3 is 
large {0{mt)), left-right mixing can also become relevant for the sbottom and 
stau. It will become clear below that At, and A,- do not contribute in a major 
way to left -right mixing, since they do not have a tan/3 enhancement. 



*Bcwarc also of different sign conventions for At- Both PYTHIA and ISAJET use the 
convention used here. 
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The Higgs boson spectrum at tree level can be expressed in terms of the 
weak gauge boson masses, the CP-odd Higgs boson mass Ma and tan/3: 1 



Ml 



hM 



1/2 M\ + Ml T -sj {Ml + M|)2 - 4M|M2 cos^ 2/3 
Ml + M^. 



(9) 



These relations yield Mh ~ Mz, but this result is strongly_modified by ra- 
diative corrections that depend on other MSSM parametersE3 The dominant 
radiative corrections to Alh grow as mf and are logarithmically dependent on 
the third-generation squark masses. The heavy CP-even and charged Higgs 
boson masses, M^ and M^±, respectively, are directly controlled by M^. If 
all SUSY particles were heavy, but M^ were small, then the low-energy the- 
ory would look like a two-Higgs-doublet model. For sufficiently large Ma, the 
heavy Higgs doublet decouples, and the effective low-energy theory has only 
one light Higgs doublet with SM-like couplings to gauge bosons and fermions. 

Within the MSSM, a general upper bound on Mh can be deter- 
mined by axareful evaluation of the one-loop and dominant two-loop radiative 
correctionsuJ For mt — 175 GeV and an extremely conservative set of assump- 
tions, f] the upper bound on the lightest Higgs mass is maximized, yielding 
Mfi ~ 130 GeV. For more moderate values of the MSSM parameters, the upper 
bound on Mh becomes smaller. Most importantly, given the general upper 
bound on Mh of about 130 GeV, the upgraded Tevatron has the potential tn 
provide a crucial test of thp minimal supersymmetric extension of the SMEj'Ej 

R-parity, defined asE^ R=(— 1)-^'^+'^^+-^, is a discrete multiplicative sym- 
metry where S is the particle spin, B is the baryon number, and L is the lep- 
ton number. All SM particles have R=l, while all superpartners have R=— 1, 
so a single SUSY particle cannot decay into just SM particles if R-parity is 
conserved. In this case, the lightest superpartner (LSP) is absolutely stable. 
Astrophysical considerations imply that a stable LSP should be color- and 
charge-neutral. The best candidates, then, are the lightest neutralino Xi ^^'^ 
the sneutrino i>, or alternatively (see below) the gravitino G. Since the LSP can 
carry away energy without interacting in a detector, the apparent violation of 
momentum conservation is an important part of SUSY phenomenologyEj Also, 
when R-parity is conserved, superpartners must be produced in pairs from a 
SM initial state. The breaking of the R-parity symmetry would result in lepton 

^In the MSSM, Ma is related to the values of the Higgs soft SUSY-breaking parameters 
and rriH^ and the Higgsino mass parameter ^ through = + "'If 2 ^ 

''To produce this bound, the masses of all SUSY particles and Ma are chosen to be around 
a TeV, tan f3 > 20, and the stop mixing parameters are varied to give the largest possible 
effect. 
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and/or baryon number violating processes. While there are strong experimen- 
tal constraints on some classes of R-parity violating interactions, others are 
hardly constrained at all. Unless it is explicitly stated otherwise, R~parity 
conservation is assumed below. 

Quite generally, the dependence of the SUSY spectrum on tan/3 can be 
very strong, and it is necessary to determine the possible range of values for 
this essential parameter of the theory. The fermion masses, which are not 
fixed by SUSY, are a function of tan/3 and the SM Yukawa couplings. For the 
up- and down-quark and lepton masses, it follows that m„ = /i„usin/3, md = 
/idWCOs/3, and mg — /i^wcos/3, where hf=u,d,e is the corresponding Yukawa 
coupling and v— 174 GeV is the order parameter of EWSB. Equivalently, 
"rriu = huV2,md = hdVi and = hivi, where V2 = wsin/3 and vi = vcosf3. 
The value of tan/3 and the Yukawa couplings can vary in a range consistent 
with the experimental values of the fermion masses. However, for the theory 
to remain perturbatively well defined up to a given cutoff scale, the Yukawa 
couplings should remain finite up to this cutoff scale. Whether this is the 
case can be determined by studying the renormalization group evolution of 
each Yukawa coupling from low to high energy scales. In particular, the large 
value of the top quark mass is associated with a large value of the top Yukawa 
coupling at low energies, which, depending on tan/3, may become^ tapJ-arge 
to be compatible with a perturbative description of the theoryCJCJ'LlEJ The 
measured value of the top quark mass, mt ~ 175 GeV, defines a lower bound 
on tan /3 of about 1.2, provided that the top Yukawa coupling remains finite 
up to a scale of the order of 10^^ GeV. If, instead, the top Yukawa coupling 
should remain finite only up to scales of order of a TeV, values of tan (3 as low 
as .5 are still possible. j] A similar situation occurs when tan/3 is large, but now 
the crucial role is played by the bottom Yukawa coupling. If tan (3 becomes 
too large, large values of the bottom Yukawa coupling are necessary to obtain 
values of the bottom mass compatible with experiment. The exact bound on 
tan /3 depends on the SUSY spectrum, since there are radiative corrections 
to the bottom mass coming from sparticle exchange loopsE3l23 Generically, it 
can be shown that values of tan/3 > 60 are difficult to obtain if the MSSM is 
expected to remain a valid theory up to scales of order 10^^ GeV. 

2. 2 Supergravity 

At present, the exact mechanism of SUSY breaking is unknown. Supergravity 
(SUGRA) models assume the existence of extra superfields (the so-called "hid- 

'This implies that a perturbative description of the MSSM would only be valid up to the 
weak scale, which is, of course, not a very interesting possibility. 
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den sector") which couple to the MSSM particles only through gravitational- 
like interactions. When SUSY is spontaneously or dynamically broken in the 
hidden sector, some of the components of the hidden sector acquire vacuum 
expectation values. Interaction terms between those components of the hidden 
sector and the MSSM superfields give rise to the effective soft SUSY-breaking 
terms of the MSSM, which are proportional vacuum expectation values of the 
hidden sector divided by powers of Afpianck- The low-energy Lagrangian then 
looks like a SUSY-conserving MSSM with a number of extra terms that break 
SUSY. Although low-energy gravitational interactions depend only on mass 
and spin, the general supergravity Lagrangian may contain higher-dimensional 
interactions between the hidden sector and MSSM superfields that are flavor 
dependent. 

The number of SUSY breaking terms is over one hundred. In the min- 
imal SUGRA scenario, however, the MSSM sparticles couple universally to 
the hidden sector, and the number of terms is greatly reduced. Using this 
guiding principle, at a scale of order Mpianck (or, approximately, Mqut, the 
scale where the gauge couplings unify), all scalars (Higgs bosons, sleptons, and 
squarks) are assumed to have a common squared-mass rriQ, all gauginos (Bino, 
Wino, and gluino) have a common mass 1^1/2, and all trilinear couplings have 
the value Aq. After specifying tan/3, all that remains is to relate the values 
of the soft SUSY-breaking parameters specified at Mqut to their values at 
Mew- This is accomplished using renormalization group equations (RGE's). 
The physical sparticle masses are then determined using relations like Eq. (^, 
or by diagonalizing mass matrices like those in Eqs. (§) and (0). Finally, 
the Higgsino mass parameter is determined (up to a sign) by demanding the 
correct radiative EWSB. At the tree level, this requires 



+ m% + Ml/2 



where tuhi and to/Zj are soft SUSY-breaking mass parameters for the two 
Higgs doublets evaluated at Mew- 

Not surprisingly, the masses of the gluinos, charginos and neutralinos are 
strongly correlated. The gaugino mass parameters Mi at the electroweak scale 
depend mainly on the running of the gauge couplings between Mqut and 
Mew- 

M3 = mi/2 ~ 2.67711/2 

acuT 

a2{M2) 

M2 = TOi/2 — 0.8mi/2 

acuT 
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Ml = "'^ ~ 0.4mi/2, (11) 

(where we consider mi/2 ~ Mew)- As wiU be shown below, once the RGE 
evolution of the Higgs mass parameters is included in Eq. (|lO|) , it follows that 
the Higgsino mass term /i tends to be larger than the Bino and Wino masses 
Mi,M2, becoming the largest for values of tan/? closest to 1. As a result, 
the lightest two neutralinos and the lightest chargino tend to be gaugino-like. 
This is similar to the example presented in Eq. (^) with the approximate mass 
hierarchy: 

M~o ~ 2M~o ~ M~± ~ l/3Mg. (12) 

A2 Ai 

The lightest neutralino x? can be the LSP. 

Because sleptons have only EW quantum numbers and the lepton Yukawa 
couplings are small, the slepton mass parameters do not evolve much from 
Mqut to Mew- The left- and right-handed soft SUSY-breaking parameters 
at the scale Mew which determine the mass eigenstates through Eq. (||) are 
given by 

"^11,2 - "^i:, - ™o + 0-5"^V2; ml^^ ~ m|;3 ~ m2 + 0.15m2/2. (13) 

For tan /3 > 40, the third generation mass parameters also receive non-negligible 
contributions from the r Yukawa coupling in their running which can modify 
these expressions. The different coefficients of mi/2 in Eq. ( p^ arise from the 
different EW quantum numbers for sleptons in SU {2)l doublets and singlets. 
If mo and toi/2 are of the same order of magnitude, physical slepton masses 
are dominated by mg. When toq is small, the sneutrino can be the LSP instead 
of the Xi- The D mass is fixed by a sum rule m?^ = mj + M^ cos 2/3. 

The squark mass parameters evolve mainly through the strong coupling to 
the gluino, so their dependence on the common gaugino mass is stronger than 
for sleptons. For the first and second generation, the left- and right-handed 
soft SUSY-breaking parameters at Mew ^re given approximately by 

"^Qi,2 - "^0 + 6.3m^/2; m^^^ ~ m|,^^ ~ml + b.Sml/^. (14) 

In general, the squarks are heavier than the sleptons or the lightest neutralino 
and chargino. Since first and second generation squark soft SUSY-breaking 
parameters are the same for squarks with the same quantum numbers, FCNC's 
are suppressed. 

For the third-generation squarks, the large top and bottom Yukawa cou- 
plings play a crucial role in the RGE evolution. As mentioned generically in 
section 2.1, the top Yukawa coupling ht is related to the top quark mass by 
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nit = (246/\/2)/if sin/3 GeV and the bottom Yukawa coupling hf, is given by 
nib — (246/-\/2)ft.fc cos /3 GeV, so that ht, is large (of order ht) when tan/3 is 
about 40 or larger. When ht at Mqut is sufficiently large, it turns out that 
its low-energy value is independent of its exact value at Mqut- This^ 
is known as the infrared fixed point solution of the top quark massE 
With the definition Yt = /ij/(47r), the infrared fixed point value of Yt at the 
scale mt is F/'' ~ Saa/O. Within the one-loop approximation, the effects of 
the top Yukawa coupling on the RGE evolution can be parameterized in terms 
of the ratio YtjY" . For small and moderate values of tan/3, the left- and 
right-handed soft SUSY-breakingjiaj|-|aineters which determine the stop and 
sbottom masses are then given byMEaEj 
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+ 





"^^3 - "^0 ( 1 - 2^ j + "^1/2 I 6-3 

^£/3 - "'0 \}-^)+ (s-S + ^ ( - V + 2^ ) ) , (15) 

and m ^.^ ~ to_Di 2 ■ For large tan /?, assuming t~b Yukawa coupling unification 
at high energies {Yi, — Yt at Mqut, which is a generic prediction of 5*0(10) 
GUT models), Xhe expressions for the third generation soft SUSY-breaking 
parameters are:l 



i6 



~ ( 1 - ^ A j + ml/^ ( 6.3 - 
2 2 2 / 1 6 Ft \ 2 



TO^^ ~TO^. ~TO^ 1- +"^!/2 5.8- -4+--^ . (16) 





-4H ^ 
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-4+l^A^ 
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Contributions proportional to Aq and A0TO1/2 with a prefactor proportional 
to (1 - Yt/Yt''^) are also present in Eqs. (|l|) and (|l|). For rrif ~ 175 GeV, 
the value of the ratio Ft/F/'" varies from 3/4 to 1 depending on tan/3, with 
Yt/Y;'- ^ 1 as tan^ ^ 1, and Ft/F/'^ ~ 0.85 for tan/? = 40. The value of At 
is governed by mi/21 and, for large values of the top Yukawa coupling, depends 
weakly on its initial value and tan /3pl 

At^(l- ^) ^0 - 2TO1/2- (17) 

The exact values of Af, and Aj- are not important, since the mixing in the stau 
and sbottom sectors is governed by the terms rub^J, tan /3 and nirH tan /3, respec- 
tively. In SUGRA models, the above relations between the mass parameters 
leads to the general prediction, to- > 0.85Mg (for the five lightest squarks and 

small or moderate tan/3). 
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The soft-SUSY breaking parameters iiiJJiii.Higgs sector also have simple 
expressions. For small and moderate tan/3, E§njC3 



Hi - '"-0 "■^"'-1/2 

k ^ ^l(^-l^) +^1/2(0.5 + ^ 1^-7 + 3^]]. (18) 



Substituting these relations back into Eq. ( pl)| ) yields the result: 



7 j tan2/3- 1 

-?/2 (0.5- (0.5+2l (-7 + 3|,)) tan^;3) (19) 

Note in Eq. (^8|) that m|^^ < 0, which is usually a sufficient condition to induce 
EWSB. For large tan/?, the Higgs mass parameters am more complicated. In 
the limit of t — 6 Yukawa unification, they simplify toEj 

- ^ ™g (1 - ^ A) + m?/, (0.5 + |f (-6 + y |,) ) , (20) 

and Eq. ( |l9[ ) must be modified accordingly. All of these relations are only 
approximate: the coefficients of toi/2 depend on the exact values of ctg and the 
scale of the sparticle masses; the coefficients of mo and Aq depend mainly on 
tan/3. 

The SUGRA model presented here is minimal (mSUGRA) in the sense that 
it is defined in terms of only five parameters at a high scale: mo, ^Oi 
tan/3, and the sign of /i. It is-natural to question exact universality of the soft 
SUSY-breaking parameters.^ For example, in a S'C/(5) SUSY GUT model, the 
left-handed sleptons and right-handed down-type squarks reside in the same 
5-multiplet of SU{5), and naturally have the common mass parameter mp^^ 
at the GUT scale. Similarly, UL,dL, ur, and cr, which reside in the same 

10-multiplet, have a common mass m'^°\ The two Higgs bosons doublets 

— ^ ^ (5') (5') 

reside in different 5- and 5-multiplets, with masses ttiq and mg . There is 

no symmetry principle that demands that all these mass parameters should 

be the same. The most naive breakdown of exact universality is to consider 

f5') (5') 

different values for uiq and itiq , taking mo as the common mass for sleptons 
and squarks. 

Depending on the exact mechanism of SUSY breaking, it may occur that 
mi/2 — 0. Low-energy gaugino masses are then dominated by contributions of 
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stop-top and Higgs-Higgsino loopsE3 In this case the gluino could be the LSP 
with a mass of Mg ~ a few GeV and the lightest neutcalino may be somewhat 
heavier due to contributions from electraweak loopsHj Light gluino scenarios 
are being explored by many experimentsa 



2.3 Gauge-Mediated Supersymmetry Breaking 

In SUGRA, gravitational interactions generate the soft SUSY-breaking terms 
in the MSSM Lagrangian. Alternatively, soft SUSY-breaking terms can be 
generated through gauge interactions. This has the feature that mass degen- 
eracies between sfermions with the same quantum numbers (and, hence, the 
same gauge couplings) occur naturally, which suppresses FCNC's. Also, in 
gauge-mediated models, the scale of the SUSY breaking is much smaller than 
the scale where gravity becomes relevant, so there is no possibility of Planck- 
scale corrections to these jiegeneracies (as them can be between the GUT and 
Planck scales in SUGRA) E3 In simple modelspl the existence of heavy messen- 
ger superfields ip with SM quantum numbers is postulated. SUSY is broken 
in a hidden sector which also couples to the messengers, so that the ip fermion 
components have mass M, while the scalar components have masses My/1 ± x, 
where x is a dimensionless parameter that controls the size of SUSY breaking. 
The MSSM gauginos and sfermions acquire masses different from their SM 
partners because of the radiative effects generated by the messenger fields. It 
is more convenient to define A = xM, which fixes A as the overall mass scale 
of the MSSM sparticles. The gaugino masses at a low energy scale /x are 

M.i^,)^^gixm~b';)A, (21) 

where i specifies the gauge group, bi is the MSSM coefficient of the beta func- 
tion for the running of a^, b'^ includes the additional effect of the messenger 
fields in the running, and g{x) ~ 1 -I- x^/6. The mass-squared of the MSSM 
scalars acquire values 



^ ^ 2 

Itt 

gixr ((6. - Kf/b,) {aUf^)/aUM) - l)] , (22) 
where Ci = 5/31^^, C2 =3/4, and C3 = 4/3 (F is the weak hypercharge) and 



f{x) ~ 1 -I- x^/36. The mass formula in Eq. (22) ignores Yukawa couplings, 
and will be modified for the stop and possibly sbottom and stau. By compar- 
ing the previous two equations, it is clear that the gaugino and scalar masses 
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are roughly of the same order of magnitude. Even after evolving these mass 
parameters to Mew (ignoring the effects of Yukawa couplings), sfermions with 
the same quantum numbers acquire the same masses, yielding a natural mass 
hierarchy between weakly and strongly interacting sfermions; the mass hierar- 
chy of the gauginos is fixed by the gauge couplings (as in SUGRA models) . If 
the superfields %p reside in a complete representation of SU{5) or 50(10), then 
unification of the gauge couplings at a high scale is not compromised, though 
the unification will occur at stronger values of the couplings and at a slightly 
different scale from the naive GUT scale. One distinctive feature of these mod- 
els is that the spin-3/2 superpartner of the graviton, the gravitino G, can play 
a crucial role in the phenomenology. Since the gravitino mass is given by the 
relation M~ — MA/Mpianck, the gravitino can be very light depending on the 
value of M, unlike SUGRA, where the gravitino has a mass on the order of 
Mw- As a result, in gauge-mediated Supersymmetry breaking, the gravitino 
can be the LSP. 

In the above discussion, it is assumed that the messengers ip form a com- 
plete GUT multiplet. However, if the messengers were neutral under some 
gauge group, then the associated gauginos would be massless at one-loop be- 
cause of gauge invariance. In particular, it is possible to construct^ model 
where the gluino is a stable LSP with a mass of a few tens of GeVEZi In this 
case, the missing energy signal for SUSY disappears, since a stable LSP gluino 
will form stable hadrons. 

2.4 R-Parity Violation 

One simple extension of the MSSM is to break the multiplicative R-parity 
symmetry. Presently, neither experiment nor any theoretical argument demand 
R-parity conservation, so it is natural to consider the most general case of R- 
parity breaking. It is convenient to introduce a function of superfields called the 
superpotential, from which the Feynman rules for R-parity violating processes 
can be derived. The R-parity violating (RPV) terms which can contribute to 
the superpotential are^ 

Wbpv = KjkL'U& + \^kVQW^ + \-^,WDW'' (23) 

where i,j,k are generation indices (1,2,3), L\ = i/j^, = and Q\ — u\, 
Q2 = are lepton and quark components of SU{2)l doublet superfields, and 
E'^ ~ ej^, — d^j^ and C/' = are lepton, down and up- quark SU{2)l 
singlet superfields, respectively. The unwritten S\J{2)l and SU(3)(7 indices 

"'In Eg. (|23[) bilinear terms are ignored. Ar^scussion of the phenomenological implications 
of such terms can be found in the literatureEJ 
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imply that the first term is antisymmetric under i ^ j, and the third term 
is antisymmetric under j <-> k. Therefore, i 7^ j in and j ^ k va 

WD^D'^. The coefficients Xijk, X^ji^ and A-^j, are Yukawa couphngs, and there 
is no a priori generic prediction for their values. In principle, Wftpv contains 
45 extra parameters over the R-parity-conserving MSSM case. 

Expanding Eq. ( ]23| ) as a function of the superfield components, the inter- 
action Lagrangian derived from the first term is 

Clle = Kjk {j>l44, + eli^le'k + + h.c.} (24) 

and from the second term, 

Clqd = {'>I44 - 6144 + 4^14 - 444+ 

(4)*44 - (4)*44 + h.c] (25) 

Both of these sets of interactions violate lepton number. The UDD term, 
instead, violates baryon number. In principle, all types of R-parity violating 
terms may co-exist, but this can lead to a proton with a lifetime shorter than 
the present experimental limits. The simplest way to avoid this is to allow 
only operators which conserve baryon-number but violate lepton-number or 
vice versa. 

There are several effects on the SUSY phenomenology due to these new 
couplings: (1) lepton or baryon number violating processes are allowed, includ- 
ing the production of single sparticles (instead of pair production), (2) the LSP 
is no longer stable, but can decay to SM particles within a collider detector, 
and (3) because it is unstable, the LSP need not be the neutralino or sneutrino, 
but can be charged or colored. 

Present data are in remarkable agreement with the SM predictions, and 
very strong bounds on the R-parity-breaking operators can be derived from the 
following processes: (a) charged-current universality, (b) T(t —^ evv)/T{T — > 
^Jil'v)^ (c) the bound on the mass of Ue, [d) neutrino-less double-beta decay, 
(e) atomic parity violation, (/) - D° mixing, (g) Re = rhadiZ°)/Te{Z°), 
(h) T{tt e9)/T{'K /zP), (i) BR{D+ n+v^)/ BR{D+ K^*e+v^), 

(j) deep-inelastic scattering, (fc) BR{t — > ttz^^-), [l) heavy nucleon decay, 
and (m) n — n oscillatioDS. Additional limits can be derived from deep inelastic 
experiments at HERAEj On the other hand, within the allowed values of the 
R-parity-violating couplings, Xijk, A^^^,, a whole new world opens up for SUSY 
searches. 



16 



2.5 Run la Parameter Sets (RIPS) 



Some CDF and D0 SUSY searches are analyzed in the framework of so-called 
"SUGRA-inspired models." To understand the limits that appear in many 
published analyses, it is necessary to state explicitly the framework behind 
RIPS. 

First, there are five main input parameters: Afg, to~, M^, tan/3 and the 
magnitude and sign of /i. The gluino mass Mg is defined to be M3, which 
is equivalent to specifying TOi/2 and, hence. Mi and M2 using the unification 



relations Eq. (11). The chargino and neutralino properties are then fixed by 
Mi,M2,tan/3 and /i. In practice, the value of /i is set much larger than Mi 
and M2, so the properties of the neutralinos, charginos, and gluino are similar 
to those in a pure SUGRA model. 

Next, all squark soft SUSY-breaking mass parameters are set to m~, and 
the Z?-terms are neglected. The result is that the first 5 squarks are degenerate 
in mass. This may be unrealistic if tan/3 is large, since the sbottom mass 
can be naturally lighter because of non-negligible off-diagonal elements in the 
sbottom mass matrix. The stop squarks are made heavier than the other 
squarks by fixing At = /i/tan/3 (see Eq. (^), which tunes away the mixing 
between and Ir- The resulting stop masses are m^^ — nii^ — y^m- + . 

Therefore, experimental limits placed on RIPS show no sensitivity to the stop 
squarks. Note that, in SUGRA models, the stop squared-mass soft SUSY- 
breaking parameters rriQ^ and m^^ are generally not equal and are smaller 
than the other squark parameters at AIew, so that one stop squark is lighter 
than the other squarks. 

Giving the other five squarks a common value at the weak scale ignores 
the details of running from the GUT scale (see Eq. (|lj)) and the different D- 
terms. However, using an average of the two formulae in Eq. (|l^, a specific 
mi/2 and m~ roughly determine a value of ttiq. Whenever there is a solution 

with niQ > (which implies m~ > 0.85A/g), RIPS has many features of a 
SUGRA model. Indeed, when m~ > Mg, the approximate SUGRA relations 
m? =m2.-0.73M|-0.27M|cos2/3, mf = m?- - 0.78A/| + 0.23MI cos 2/3, 
and m? — m~ — 0.73Af| + 0.5M|cos2/3 are used to fix the slepton masses.^ 
The region m~ < A4g is very hard to realize in SUGRA models, but is also 
worth investigating. In this case, for some analyses, a constant value of 350 



"Observe that the D— terms for the sleptons, although correct, are negligible in comparison 

with the approximation made in defining a common m~. However, D-terms are included 

_ Q 

to assure the correct splittings between the II and v masses. 
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GeV is set by hand for nij^, m^^, and mp. Accordingly, experimental limits 
placed on RIPS when m~ < Mg show little sensitivity to the sleptons. 

Finally, the Higgs mass Mj^ is used to determine the Higgs boson sector. 
This is equivalent to considering partial non-universality for the scalar sfermion 
and Higgs boson soft SUSY-breaking mass parameters at high energies, i.e. 
Too 7^ TOq'^ ^ TOq^ (see the discussion near the end of Sec. 2.2). In practice, 
the CP-odd Higgs boson mass Ma is set to a large value, so that the lightest 
neutral Higgs boson h has SM-like couplings to gauge bosons and fermions, 
and all other Higgs bosons are so heavy they are not kineniatically accessible 
at the Tevatron. 



3 The CDF and D0 Detectors 



The CDfEI and 000 detectors are located at the interaction regions BO and 
DO in the accelerator ring:| Both detectors feature particle tracking detectors 
close to the interaction region, surrounded by quasi-hermetic calorimetry cov- 
ering the region of pseudorapidity^ of approximately \ri\ < 4. Muon detection 
systems are located outside the calorimeters for both detectors. 



3.1 The CDF Detector 

The CDF detector is distinguished by its magnetic spectrometer: a 3-m di- 
ameter, 5-m long superconducting solenoidal magnet, which creates a 1.4 T 
field uniform at the 0.1% level and contains the parlicle tracking detectors. A 
four-layer silicon microstrip vertex detector (SVX)p3 located directly outside 
the beampipe, tracks charged particles in the r — </> plane. The SVX measures 
the impact parameter of tracks coming from secondary vertices of bottom and 
charm decays with a typical resolution of 30 /um, providing heavy-flavor tag- 
ging for jets. A set of vertex time projection chambers (VTX) surrounding the 
SVX provides tracking in the radial and z directions and is used to find the z 
position of the pp interaction. Outside the VTX, between from a radius of 30 
to 150 cm, the 3.2-m long central tracking chamber (CTC) is used to measure 
the momentum of charged particles, with up to 84 measurements per track. 

The calorimeter is divided into a central barrel (|?7| < 1.1), end-plugs'! 
(1.1 < \r]\ < 2.4) and forward/backward modules (2.4 < \ri\ < 4.2). Each of 

"The Tevatron ring has six-fold symmetry, with the centers of the straight sections la- 
belled as AO, BO, CO, DO, EO, and FO. 

''The pseudorapidity rj is defined as — ln(tan ^). In the CDF and D0 coordinate systems, 
9 and <p are the polar and azimuthal angles, respectively, with respect to the proton beam 
direction z. 

"End-plugs" because they plug into the ends of the solenoid and central calorimeter. 
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these is segmented into projective ^ electromagnetic and hadronic towers sub- 
tending 0.1 in 77 by 15° in (j) in the central calorimeter and 5° elsewhere. 
Wire chambers with cathode strip readout give information on electromagnetic 
shower profiles in the central and plug calorimeters (CES and PES systems, 
respectively). A system of drift chambers (CPR) outside the magnet coil and 
in front of the electromagnetic calorimeters serves as a "preradiator," allowing 
additional photon/7r° discrimination on a statistical basis. Muons are identi- 
fied with the central muon chambers, situated outside the calorimeters in the 
region \r]\ < 1.1. 

The magnetic spectrometer measures muon and other charged particle 
transverse momenta with a resolution CTp^/PT < O.OOlpr (pt in GeV) and al- 
lows a precision calibration of the electromagnetic calorimeters by comparing 
the measured calorimeter response to the measured momentum from high- 
energy electrons from W decayscJ Electron energies are measured with a reso- 
lution aE/E = .135/v^® .01 (Et in GeV).{] Jets arareconstructed as energy 
clusters in the_calorimeter, -using a cone algorithm □ with a cone radius of 
either R = 0.7^ or i? = 0.4EI in Ar] x A(j> space. The jet energy resolution in 
the central region is approximately as/E = .80/^/E (B .04 {E in GeV). 

Missing transverse energy,[] (.^y), a key quantity in SUSY searches, is 
calculated as ^ £;tower . ^ where the sum is over both electromagnetic 
and hadronic calorimeter towers^ in I77I < 3.6, i<;tower |-j^g energy measured 
in the tower, n is the unit vector pointing in the direction of the center of the 
tower from the event vertex, and f is the unit vector in the radial direction. 
The l^rp is always corrected for the momentum of muons; for many SUSY 
analyses, it is also corrected for the calorimeter response to jets. The typical 
resolution on a component of is 5.7 GeV in Z'^ — s- e+e^ events. 

3.2 The D0 Detector 

The D0 detector consists of three major components: a non-magnetic central 
tracking system, central and forward Liquid Argon sampling calorimeters, and 
a toroidal muon spectrometer. The central tracking system consists of four 
detector subsystems: a vertex drift chamber, a transition radiation detector, 
a central drift chamber, and two forward drift chambers. Its outer radius is 

'^Projective means pointing approximately at the interaction region. 

■^The symbol © denotes addition in quadrature, e.g. aS) b = \/a^ + b'^. The total resolu- 
tion can be parameterized this way when there are two or more independent components of 
resolution. 

•^The transverse momentum of a particle with momentum p is pT = p sin 9. The analogous 
quantity using energy, defined as Ex = Esm9, is called transverse energy, 
tower is a cell in jj — <^ space. 
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76 cm. The system provides identification of charged tracks in the pseudora- 
pidity range \r]\ < 3.5. It measures the trajectories of charged particles with a 
resolution of 2.5 mrad in (jj and 28 mrad in 9. Using the reconstructed charged 
tracks, the position of the primary interaction along the beamline direction 
is reconstructed with a resolution of 8 mm. The central tracking system also 
measures the ionization of tracks to allow single charged particles to be dis- 
tinguished from e+e^ pairs from photon conversion. The transition radiation 
detector aids in distinguishing electrons from charged pions. 

The calorimeter is transversely segmented into pseudoprojective towers 
with A?7 X A(j) — 0.1 x 0.1 and provides full coverage to I77I < 4.2. The 
calorimeter is divided into three parts, a central calorimeter and two end 
calorimeters. These are further segmented into an inner electromagnetic sec- 
tion, followed by a fine hadronic section, and then a coarse hadronic section. 
Between the central and end-cap calorimeters, a set of scintillator tiles pro- 
vides improved energy resolution for jets that straddle the two detectors. The 
electromagnetic (EM) calorimeter is divided into 32 modules in (p, each of 
which has 22 layers, each approximately 1 radiation length ^ thick, with Liq- 
uid Argon as the active element and 238[/ plates as the passive element. These 
layers are arranged into four longitudinal segments per tower, called cells. 
The first cell contains 2 layers, the second cell contains 2 more layers, the 
third cell is finely segmented, with Ary x Ai^ = 0.05 x 0.05 and contains 
7 layers, and the last cell contains 10 layers. The fine hadronic calorime- 
ter uses a Uranium-Niobium alloy as its passive element, and the coarse 
hadronic uses copper. The electron energy resolutions, as measured in the 
EM calorimeter, are as/E = 0.130/\/;B^ © 0.0115 © 0.4/K for I77I < 1.1, and 
aE/E = 0.157/v^© 0.010© 0.4/i; for 1.4 < |?7| < 3.0H The azimuthal po- 
sition resolution for electrons above 50 GeV as measured by the calorimeter 
is 2.5 mm. The muon spectrometer provides muon detection in the range 
|?7| < 3.3. The total thickness of the calorimeter plus the toroid varies from 13 
to 19 interaction lengths, making hadronic punch-through backgrounds negli- 
gible. The muon momentum resolution is cTp/p = 0.18(p — po)/?' © O.OOSp {p 
in GeV/c, po = "2 GeV/c). The D0 detector has a compact tracking volume 
which helps control backgrounds to prompt muons from in-flight decays of tt 
and K mesons. Jets are reconstructed as energy clusters in the calorimetet, 
using a cone algorithm with a cone radius of i? = 0.5 or i? = 0.7 in Ary x A^^ea 
The jet energy resolution is approximately oe/E = 0.8/^/E {E in GeV)E2l 

]^rp is calculated using the vector sum of energy deposited in all calorimeter 
cells, over the full calorimeter coverage for \r]\ < 4.2, with corrections applied 

''The radiation length is the mean distance traversed by an electron in a given material 
during which it radiates a fraction 1 — of its energy via bremsstrahlung. 
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to clustered cells to take account of the jet energy scale, and to unclustered 
cells as determined from studies of Et balance in e+e^ events that do 

not contain hadronic calorimeter clusters. The resolution on a component of 
]^rp in "minimum-bias" events] is 1.1 GeV+ 0.02 J2 ^t, where ^ Et is the 
scalar sum of transverse energies in all calorimeter cells. For some analyses, 
the $!rp is corrected for the presence of muons, which only leave a small fraction 
of their energy in the calorimeter. 

3.3 Experimental Realities 

There are potentially two types of backgrounds to any experimental signature, 
physics and instrumental. Physics backgrounds mimic the event signature 
even in an ideal detector, while instrumental backgrounds arise because of 
detector flaws. Experimental signatures - SUSY or otherwise - are identified 
from "objects" - the building blocks of the event. Examples of objects that 
CDF and/or D0 use are electrons, muons, tau's, fj,, so-called "generic jets" 
(presumably from quarks and/or gluons, but without flavor identification), c's, 
&'s, photons, and, using another level of kinematic reconstruction, W^'s, Z's, 
and Vs. The selection of each of these objects carries with it an efficiency and 
also a "fake rate," a probability that the object is actually a different object 
which has been misidentified. The description of an object as an "electron," for 
example, more precisely means an "electron candidate that passes the electron 
cuts," and no more and no less. 

For the majority of searches, the signal and its physics backgrounds can 
be estimated using Monte CacIo simulation. The output of an event generator 
such as ISAJETEI or PYTHIAB can be folded with relatively simple parameteri- 
zations of the detector response to give a good description of the data. A typical 
simple simulation transforms the final state partons from a Monte Carlo into 
jets, using a clustering algorithm similar to the one used for the data. It then 
convolutes the momenta of the electrons, photons, muons, and jets with the 
appropriate experimental resolutions, generating "smeared" momenta. is 
calculated by first summing the smeared visible momenta, and then adding the 
effects of additional minimum-bias events in the same beam crossing. When 
calculating the geometric acceptance of the detector, it is necessary to include 
the distribution for the interaction vertex position in z, which is Gaussian with 
an RMS of approximately 30 cm. Finally, the detection efficiencies for elec- 
trons, photons, 6-quarks, r's, and muons (jets above 20 GeV are usually found 

'These events are defined by the requirement that a beam-beam collision took place as 
measured by arrays of scintillation counters forward and backward near the beampipe, and 
therefore have a smaller selection bias than events selected with more selective triggers. 
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with good efficiency) are applied. Initial and final state gluon radiation need 
to be included, since they can affect the efficiency by adding extra jets which 
can modify the event signature or "promote" backgrounds into the signature. 

In order to make even rough predictions of instrumental backgrounds, im- 
perfections in the detectors must be taken into account. Two effects make these 
difficult to estimate: fakes, and tails on jet energy resolution distributions. Be- 
cause of the very large multi-jet production rate at the Tevatron, there can be 
significant fake backgrounds, even if the fake probability is very small. Fakes 
are very complicated, and the fake rate must be evaluated for each analysis 
using the appropriate data. In general, the efficiency for properly identifying 
an object and the probability that another object fakes it are complementary. 
For signatures dominated by instrumental backgrounds, tighter selection crite- 
ria ("cuts") make for a purer sample, but reduce the efficiency. For signatures 
dominated by physics backgrounds, looser cuts are preferred because they pro- 
duce a higher efficiency for the same ratio of signal to background. Fakes are 
an especially serious problem for signatures involving photons, tau's, b- and c- 
quark tagging, and ]^rp in events with jets. Although jet energy resolutions are 
roughly Gaussian, even small non-Gaussian tails, convoluted with the large jet 
cross section, can lead to significant numbers of events with large fake ]^rp. In 
addition, there are some other factors which contribute to fakes and which are 
unique to working at the Tevatron, such as the long interaction region, the ex- 
istence of multiple collisions in a single event, the presence of the Main Ring in 
the same tunnel as the beam, and larger cosmic ray backgrounds than found at 
detectors that are deeper underground (such as the LEP experiments). Even a 
full detector simulation cannot correctly model all detector imperfections and 
these other effects. 

4 The Present Status of Sparticle Searches 

Most previous Tevatron searches have been made under very specific assxixortj 
tions. Several of the classic-signatures , such as "jets-|--^2.,"E3 "trileptons,"oEj 
and "same-sign dileptons"LJ are likely to be fruitful in many models; others 
may be specific to a certain model. We advocate a signature-based approach, 
in which a broad range of channels are studied for departures from SM ex- 
pectations, without engineering the analysis for a specific class of models (see 
Appendices A and B). While this may sound obvious, it is a large task with no 
well-defined beginning. With the experience of several years of data taking, 
however, experimentalists now have an idea of what they can do well. Sig- 
natures involving high-pT isolated electrons, muons, and photons, 6-quarks, 
c-quarks, r leptons, and/or can be measured accurately and have under- 
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standable backgrounds. Such signatures then are a practical starting point 
for the new generation of searches. In addition, there are motivations that 
are more general than the predictions of specific models for studying samples 
of (i) high-mass, high-,^^ events to probe gluino and squark production, (ii) 
inclusive leptons, lepton pairs, and gauge bosons (7, W, and Z) to probe both 
direct and cascade production of charginos and neutralinos, and (in) third- 
generation fermions {t, b, and r) to probe decays of light squarks and sleptons, 
as well as decays of Higgs bosons and Higgsinos. When setting limits, it is 
convenient to use specific models which reduce the number of free parameters 
(such as SUGRA or RIPS), but the quoted limits are valid only in that context 
and are not general limits. New physics may appear where we do not expect 
it. 

In the sections below we discuss the present status of searches at the Teva- 
tron. We also discuss the phenomenology behind these searches and comment 
on possible improvements. Table |l| summarizes those CDF and D0 analyses 
that have been published or presented at conferences. The pace and scope of 
supersymmetry searches at the Tevatron, as well as the sophistication, have 
grown enormously in the last several years as the emphasis has shifted beyond 
the top quark and more data have become available; there are many analyses 
currently in progress. A much broader picture of the Tevatron's capabilities 
should emerge as these results become available. 



4-1 Charginos and Neutralinos 

In SUGRA models, the light neutralinos and charginos are much lighter than 
the gluino or squarks, and may be the only sparticles directly accessible at 
the Tevatron. In general, the lightest neutralino is a good LSP candidate, so, 
assuming all charginos and neutralinos are relatively light, a discussion of their 
phenomenology is a good starting point for an overview of Tevatron searches. 

Chargino and neutralino pairs would be produced directly^ at hadron col- 
liders through their electroweak couplings to squarks and the vector bosons 7, 
W, and Z. The production cross sections are not a simple function of chargino 
and neutralino masses, but depend also on their (model-dependent) mixings 
and the squark masses. Quite generally, there are three contributions to XX 
production: (i) s-channel gauge boson production, (ii) t-channel squark ex- 
change, and (Hi) interference. For type (i), the reactions qq — > W* xtx'j 
and qq /Z* xtxj occur through Wino and Higgsino components, 

and qq ^ Z* ^ XiX'j through Higgsino components of the neutralinos and 

"They may also be produced indirectly in the decays of heavier sparticles. 
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Table 1: A compilation of results from Run I Tevatron SUSY searches as of the summer 
of 1997. The symbol b denotes an additional fe-tagged jet. Also listed are the references 
and the section of this chapter where each analysis is discussed. More info rmation is avail- 
able for D0 at Ijittp : //www — dQ.f iial.aov I'Dubliclnew Ine^ u _puhlic.ht'm^ and for CDF at 



http : / /www — cdf.fnal.gov/ 



Sparticle 


Signature 


Expt. 


Run 




Ref. 


Sec. 


Charginos 
and 

Neutralinos 


_^2^-|-trilcpton 


CDF 


la 


19 




M 




4.1 


_^2^-|-trilcpton 


CDF 


lab 


107 




^ 1 
3i 






.^y+trilepton 


D0 


la 


12.5 




^ o 
3Z 






^7^+trilepton 


D0 


lb 


95 




30 






'7'7+-^^or jets 


CDF 


lb 


85 




1 r 

Ah 


1 


4.10 




D0 


lab 


106 


1 


or 


1 
1 




Squarks 
and 

Gluinos 


]^rp+ >3,4 jets 


CDF 


la 


19 




38 




4.2 


>3,4 jets 


D0 


la 


13.5 




37 




■)■) 


>3 jets 


D0 


lb 


79.2 




39 




■}■) 


dilepton+ >2 jets 


CDF 


la 


19 




74 




■)■) 


-^y+dilepton+ >2 jets 


CDF 


lb 


81 




75 




■)■) 


_^y+dilepton 


D0 


lb 


92.9 




73 






Stop 


t^T+^+ >2 jets+6 


CDF 


lb 


90 




30 




4.3 


^T+^+ >3 jets+6 


CDF 


lab 


110 




32 




75 


dilepton+jets 


D0 


lb 


74.5 




31 




■)■) 


^j,+2 jets 


D0 


la 


7.4 




79 








CDF 


lb 


85 




HE 


1 


4.11 


Sleptons 


77 


D0 


lab 


106 




12C 


1 


4.10 


Charged 
Higgs 


dilepton + 


CDF 


la 


19 


1 


IOC 


1 


4.5 


T+2 jets+^j. 


CDF 


la 


19 








77 




CDF 


lab 


91 




37 




77 


T+6+^y+(f,T,jet) 


D0 


lab 


125 








77 


Neutral 
Higgs 


WH -^£ + ^r+fe+jct 


CDF 


lab 


109 




35 




4.6 


WH -^£ + ^j^+fe+jet 


D0 


lb 


100 




M 




77 


WH, ZH 77+2 jets 


D0 


lb 


101.2 




37 






ZH 6+jet+^j, 


D0 


lb 


20 


|101 


1 


77 


WH, ZH ^2 jets+2 b's 


CDF 


lb 


91 


fee 






R violating 


dilepton+> 2 jets 


CDF 


lab 


105 


1 


104 


1 


4.8 


Charged LSP 


slow, long-lived particle 


CDF 


lb 


90 




102 


1 


4.7 



charginos, respectively.|| In type (ii), the scattering quarks and antiquarks 
exchange squarks subject to the constraint that Qn couples only to Bino com- 

''An asterisk superscript * on a gauge boson refers to a resonance off the mass shell. 
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poncnt of ncutralinos, and Ql to Bino and Wino components of the neutralinos 
and to charginos. Any Q/? or Ql eoupling to Higgsino-components is propor- 
tional to the corresponding quark mass and can be ignored for chargino and 
neutralino pair production (even though h and i have large couplings, the con- 
tribution of h and t initial states from hadrons is small). Hence, in chargino 
pair production {e.g. XiXi)i only Ql exchange is important, since Qu only 
couples to the charged Higgsino. In the case when all squarks are heavy, type 
(«) contributions dominate. When the squarks are sufficiently light, both type 
{a) and {in) can be important. For example, if ^ Mi,M2 (see Eq. (||)), 
then X?, X2' ^'^d xf are mostly gaugino-like, and both Ql and Qr can have 
electroweak strength couplings to them. If, on the other hand, <C Mi,M2 
(see Eq. (^), Q exchange is only important for the heavier states, which might 
not be kinematically accessible. 

Figure |^ shows the production cross sections of various chargino and neu- 
tralino pairs at the Tevatron for the limiting cases considered earlier in Eq. (^ 
and Eq. (||). In this figure, tan/3 = 2, m~ = 500 GeV, and the gauginos obey 
the unification relations (see Eq. (|ll|)). The left figure is generated by fixing 
/i = —1 TcV and varying mi/2, the right figure by fixing mi^2=^ TeV and 
varying /i. For reference, the SM W^h production cross section is shown as a 
function of Mh. The example of Eq. (||) is most like a pure SUGRA model, in 
which the couplings of WX1X2 ^^'^ ^xtxi ^-re large, so the xf X2 ^^'^ xtxi 
production cross sections are the largest. 

Charginos and neutralinos can also be produced in associated production: 
Qx and gx- This is discussed further in Section 4.2. 

The decay patterns for the charginos and ncutralinos are also very model 
dependent. When kinematically allowed, a tree-level 2-body decay domi- 
nates over a tree-level 3-body decay, because the latter has an extra factor 
of (7^/(47r^) in the decay rate. Possible 2-body decays of the chargino are to 
W^Xi^ H^Xi^ ^L.Rf, v£ and Qq' . The heavier chargino X2' can also decay to 
Zxf and hxf (when Higgs bosons are part of the event signature, the final 
states can contain heavy flavor quarks or tau leptons). When no 2-body final 
states are kinematically allowed, the chargino will decay to a 3-body final state 
with contributions similar to those for chargino production (described previ- 
ously): (i) virtual gauge boson decays, {ii) virtual sfermion decays, and {in) 
interference. A common decay is x^ ~^ X^ff'- If sfermions are much heavier 
than a chargino, then type (?) dominates, and the decays proceed through W* 
with branching ratios similar to those of the on-shell W boson. Virtual squarks 
and sleptons can significantly alter the branching ratios to a specific //' final 
state depending on the squark and slepton masses, so that a 100% branching 
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Figure 1; Production cross sections at the Tevatron for chargino and neutralino pair produc- 
tion versus the lightest chargino mass for two Umiting models discussed in the text: large 
(Left) and small |//| (Right). The Wh cross section (curve 6) is shown for reference as a 

function of Mf^. 



ratio to ifXi or jjXi final states is possible. A 3-body decay xf 9QQ' can 
occur through virtual squark decays, provided the gluino is light enough (see 
Appendix A). 

Similarly, J'or neutralinos, the 2-body decays to Zx'^, /ix", W^^X^, H^x^, 
(■l,r(-, or Qq will dominate when kinematically allowed. When this is not 
the case, the 3 body decay x'jff (or X*^ ^ SIQ) can occur through virtual 

Z bosons, squarks, or sleptons. Three-body decays can also be in competition 
with a loop decay, so that x^ — > Xil^ since the same factor of g^/ (47r)^ coming 
from a loop integral for a 2 body decay also appears in the 3 body decay rate. 
Such a decay is important when the Xi is Higgsino-like and X2 is gaugino-like, 
or vice versa. 

In general, for X2xt production, the final states are (i) four leptons and 
or (a) two leptons, two jets and l^j., or (iii) four jets and $!rp. Some of 
the leptons can be neutrinos. For xtxi production, the final states are (i) 
two acollincar charged leptons and (ii) one charged Icpton, 2 jets, and 
l^q,. or (Hi) four jets and A wide variety of signatures is possible from the 
production of other chargino and neutralino combinations. 
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Trileptons 

The production of xf X27 followed by the decays xt Xi^^ and X2 ~^ ^^^~Xij 
is a source of three charged leptons (e or /i) and called trilepton events 
(the Ifrp is silent). The trilepton signal has small SM harkgrounds, and is 
consequently one of the "golden" SUSY signatures OlZIC3'E3 

The overall efficiency for xtx2 production with decays into three detected 
leptons is set mainly by the branching ratio for the trilepton final state, which 
is highly model-dependent. The efficiency depends on mass splittings between 
the Xi S'lid. X2 ai^d the Xi- For example, if the 2-body decay chain X2 ~^ 
and £ £xi occurs, and the mass splitting between X2 ^^.d i or between £ and 
Xi is small, one of the leptons can be too soft to detect. In addition, if the 
mass splitting between Xi a-nd X2 large, decays to real Z or ft, bosons are 
possible. Real Z bosons have a small branching ratio to e'^e^ and fJ-^ fJ-~ and 
are a SM background, and h will decay mainly to bb, decreasing the trilepton 
rate. A similar discussion holds for the decays of the xt j especially in SUGRA 
since M~± ~ ^^o- For example, in SUGRA models, the branching ratios can 

be BR{xt £iyx'i) = 0.22,Bi?(x^ £~^£~Xi) = 0.32 when the sleptons are 
off-shell but lighter than the squarks, or BR{xf — > £1^X1) — 0.66, BR{x2 
£'^£~Xi) — when the sneutrinos are light enough to allow X2 ~^ on-shell. 
The branching fractions also depend on tan/3, and, for large tan/3, decays to 
Vs and r's are enhanced. The decays of X2 to e^e^ and /i^/i^ are strongly 
suppressed for large tan/3, faUing a factor of about 5 bet-ateen tan/3=2 and 
tan/3=20 if the squarks are much heavier than thesleptonsEj 

The results of the CDulio and D0 searches Mrn are shown in Fig. an- 
alyzed using RIPS (see Sec. ^.5| ). The searches include four channels: e+e~e''=, 
e+e"//^, e^/i+/i~ and /z^/x+/Lt~. The CDF analysis is based on the cuts listed 
in Table |[ and requires one lepton with Et >11 GcV, passing tight identifi- 
cation cuts, and two other leptons with Et >5 GcV (electrons) or px >4 GeV 
(muons), passing loose identification cuts. All leptons must be isolated, mean- 
ing there is little excess Et in a cone of size R = 0.4 in 77 — space centered 
on the lepton. The event must have two leptons with the same flavor and 
opposite sign. If two leptons of the same flavor and opposite charge have a 
mass consistent with the J/^, T or Z boson, the event is rejected. After this 
selection, six events remain in the data set, while the expected background, 
dominated by Drell-Yan pair production plus a fake lepton, is 8 events. After 
demanding ^7i> 15 GeV, no events remain, while 1.2 are expected from SM 
model sources. 

The D0 analysis requires leptons with Et > 5 GeV satisfying the selection 
criteria of Table 0. However, several different triggers are used, and some 
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Figure 2: (Left) The CDF 95% C.L. limits on cross section X branching ratio for X2 
production in 107 pb~^of data. The hmit is on the sum of the final states eee,eeii, nfifi 
and fifie when xj*^ ^i^Xi s-i^d X2 rH^^X?- The signals expected for three different RIPS 
scenarios are shown for comparisonEZI Typically, M~g ~ ~ 2M~q . (Right) Similar 

X2 Xi X-i 

limits from D0, but for the average of all four channels. The curves (A), (B), and (C) show 
the Run la, Run lb, and combined limits. Curve (i) shows the predicted cross section X 
branching ratio assuming BR(xJ — > ^1^X1)= BR(x2 — > ^~''^~Xi)=l/3 = e,^, r). Curve 
(ii) assumes BR(x^ ii/x1)=0.1 and BR(x2 ^ £+€-x?)=0.033. For both CDF and D0, 
kinematic efficiencies are calculated using the production cross section from ISAJET. 



lepton categories are required to have a larger Et to pass the various trigger 
thresholds. All leptons are required to be isolated. To reduce events with 
mismeasured $!rp^ the Jfj, must not be along or opposite a muon. Additional 
cuts arc tuned for each topology. For example, the background from Drell- 
Yan pair production plus a fake lepton is highest in the eee channel, so these 
events are rejected if an electron pair is back-to-back. The ]^rp cut is 15 GeV 
for eee, and 10 GeV for the other three topologies. No events are observed in 
any channel with a total of 1.26 events expected from (i) Drell-Yan production 
plus a fake lepton and (ii) heavy-flavor production. 

To compare the D0 and CDF 95% C.L. results, note that the two exper- 
iments present different quantities: the D0 limit is on the "average" of the 
4 modes (eee, ee/x, e/i/j,,and /i/i/i), while the CDF limit is on the sum. After 
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Table 2: Selection criteria and results for the CDF trilepton gaugino search. The 'Very Loose' 
muon category refers to isolated stiff tracks that leave only small amounts of energy in the 
calorimeters, but do not have a corresponding track in a muon chamber (this substantially 

increases the acceptance). 



Quantity 


Criteria/Cut Value 


Lepton ID Categories 


1 Tight + 2 Loose, jSQI ^ 3 


Lepton Isolation 


EEt < 2 GeV in a cone R = 0.4 


Tight E^, T] range 


> 11 GcV, \r]\ < 1.0 


Loose E^, T] range 


> 5 GeV, \ti\ < 2.4 


Tight p^, 1] range 


> 11 GeV, Irj] < 0.6 


Loose pt^ , rj range 


> 4 GeV, \r]\ < 1.0 


Very Loose Muon p^, rj range 


> 10 GcV, \r)\ < 1.0 


ZjTjJ/ip mass window cuts 


75- 105,9- 11,2.9- 3.1 GeV 


A(j) between highest 2 Et leptons 


< 170° 


AR between any 2 leptons 


> 0.4 




> 15 GeV 


jCdt 


107 pb-^ 


Expected Background 


1.2 events 


Observed Events 


events 



Table 3: Selection Criteria and Results for the D0 Run lb Trilepton Search. 





Channel and Trigger 




eee Trigger 


ee/i Trigger 


e/i/i Trigger 


MMM Trigger 


Energy ordered 


>22,5,5 e$T 


> 22(e), 5, 5 e$T 


> 9(c), 10, 5 efi 


> 17,5,5 ^i 


Et (GcV) 


> 14, 9, 5 2e$T 


> 14, 9, 5(Ai) 2e$j, 


>17(/i),5,5 /X 


> 5, 5, 5 /i/i 






> 9(e),10(/i),5 en 


> 5, 5, 5 /x/i 




Mass window cut 


80-100 GeV 


None 


0-5 GeV 


0-5 GeV 




> 15 GeV 


> 10 GeV 


> 10 GeV 


> 10 GeV 


Acc- 


ItT- A0e,e| >0.2 


None 


Itt - A<?i^,^| >0.1 


Itt- A(/.^,^| >0.1 


ents 


2 leading e's 






all combinations 


jCdt 


94.9 pb-^ 


94.9 pb-^ 


89.5 pb-^ 


75.3 pb-^ 


Background 


0.34 ±0.07 


0.(il±0.3G 


0.11 ± U.04 


0.20 ± 0.04 


Observed 















accounting for this difference, the CDF limit is twice as sensitive at a given 
xf mass. The CDF limit shown is compared to three RIPS, which have differ- 
ent ratios of m~ to Mg. The D0 limit is compared to a wide variation of pos- 
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sible branching ratios. Curve (i) assumes BR(x5^ ih'Xi) — BR(x2 ^ ^^Xi) 
= 1/3 (no hadronic decays), while curve (m) assumes BR{xi ^vxi) = 0.1 
and BR(x2 ^^Xi) = 0.033 (gauge boson-Uke decays). The D0 theory curve 
assumes heavy squarks, suppressing the squark exchange diagram but the CDF 
curves do not. The wide differences in the theory curves in Fig. show the 
dangers of quoting a mass hmit rather than a cross section x branching ratio 
hmit. 




40 50 70 100 200 300 400 500 

Lightest Chargino Mass (GeV) 



Figure 3: The overall trilepton signal rateO for an ensemble of SUGRA models as a function 
of the lightest chargino mass M~-j- . The kinematics cuts are different than those used in the 

^1 

present experimental analyses. The different symbols refer to solutions showing interesting 
behavior where X2 {^) ^ neutral "invisible" branching ratio (generally X2 ~* '^^) > 90%, 
(B) a large destructive interference in 3— body leptonic decays, (C) a branching ratio to Higgs 
> 50% dominates, or (D) all other solutions. The horizontal lines represent the reach for 

various integrated luminosities. 

The experimental limit on the cross section (times branching ratio) de- 
pends on the kinematics of the decays, mostly through the mass splitting be- 
tween xf, X2 LSP. For /i > 0, the mass splitting M~± —M~o is smaller, 
and the lepton pr cuts are less efficient. As long as the sleptons are heavier 
than xt and X2j the leptonic decays of xt and X2 through virtual W and Z 
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bosons and slcptons have similar kinematics since the decays are dominated 
by phase space. However, when the experimental result is presented as a limit 
on the mass of the lightest chargino rather than as a cross section limit, the 
result is highly model-dependent. The theoretical cross section and branching 
ratios are strongly affected by the SUSY parameters. If sleptons and squarks 
are both heavy, the decays of xt and X2 to leptons follow the pattern of the 
SM gauge particles. If the sleptons are light and the squarks are heavy, decays 
to trileptons are enhanced. 

Figure H shows the wide variation in total efhciency x cross section x branching 
ratio at a given chargino maas-for the trilepton signature by sampling a large 
ensemble of SUGRA modelsO The horizontal lines represent the reach for 
various integrated luminosities, showing that the Tevatron reach can be quite 
good for sufhciently high luminosities. However, even in the restrictive SUGRA 
framework, no absolute lower limit on the chargino mass is possible. 



4-2 Squarks and Gluinos 

Since the Tevatron is a hadron collider it can produce gluinos and squarks 
through their SU{3)c couplings to quarks and gluons. The dominant produc- 
tion mechanisms are gg, qq gg or QQ*, qq QQ and qg — > Qg, qg Q*g- 
Because QCD is unbroken,"! the production cross sections of gluinos and squarks 
can be calculated as a function of only the squark and gluino masses (ignoring 
EW radiative corrections). Figure ^ shows the production cross sections for 
squarks and gluinos as a function of the sparticle masses at ^/s = 1.8 TeV (left) 
and 2 TeiZ. (right), where NLO Supersymmetric QCD corrections have been 
includedE3 The total cross sections can be of the order of a few picobarns for 
squark and gluino masses up to 400 GeV. The NLO corrections are in general 
significant and positive (evaluated at the scale^Q — rh, the average mass of the 
two produced particles), and much less sensitive to the choice of scale than a 



"The strong couplings of gluinos and squarks are the same as those of gluons and quarks, 
so that the production cross sections are the usual strong interaction cross sections. 

^ In the perturbative calculation of scattering probabilities in field theory, two scales 
appear: a factorization scale, where the parton distribution functions are evaluated, and 
a renormalization scale, where the running coupling is evaluated. In practice, these scales 
are chosen to be the same. This scale should be representative of the typical momentum 
flowing through a Feynman diagram. In Drell-Yan production, for example, the scale is the 
invariant mass of the Drell-Yan pair. The higher the order of a perturbative calculation, the 
lesser the dependence on this scale. 
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LO calculation.^ In Fig. ^ five degenerate squark flavors are assumed,^ which 
is one way to suppress FCNC's. 




Figure 4: (Left) Production cross sections for gluinos and squarks versus sparticle mass Mx 

at the Tevatron, y/s = 1.8 TeV, assuming degenerate masses for 5 flavors of squarks. For 

QQ* and QQ production, Mx is the squark mass and Afg = 200 GeV. For gg production, 

Mx is the gluino mass and A/~ = 200 GcV. For Qg production (a), Mx is the squark mass 

and Mg = 200 GeV; for (6), Mx is the gluino mass and M~ = 200 GeV. The scale is the 

Q 

average mass of the two produced sparticles. (Right) The same curves with y's = 2 TeV. 
The bands show the change in rate from varying the scale from 1/2 to 2 times the average 

mass of the produced particles. 



The ultimate detectability of squarks and gluinos depends upon their de- 
cays, which, in turn, depends on the electroweak couplings of the squarks and 
the mixings in the neutralino and chargino sector. Since squarks and gluinos 
decay into charginos and neutralinos, their signatures can be similar to XX pro- 
duction, but with accompanying jets. If m— > Mg, then the squark has the 

2-body decay Q gq. The gluino has then the possible decays g qqXi or 
g qq'xfj where q can stand for t or 6 as well, or even g ti* or ti if kinemat- 
ically allowed. The gluino can also decay via one-loop diagrams as g ^ gXi- 
If, instead, m~ < Mg, then the gluino has the 2-body decay g — > Qq. The 

The inclusion of NLO effects in the cross sections will typically raise the lower bounds 
for squark and gluino masses by 10 to 30 GeV with respect to the LO cross sections evaluated 
at a scale equal to the invariant mass of the produced particles. 

Assuming the same mass for the bottom squark is a simplification that becomes ques- 
tionable in the large tan /3 region. 
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squarks can then decay as Ql,r QXi^ ul ^xt ^ ^^"^ "^xT ■ The final 

event signatures depend on the decay channels of the charginos and neutrali- 
nos, but, typically involve l^j^ and a multiplicity of jets and/or leptons. The 
3-body or top-stop decay modes of the gluino can produce a higher multiplicity 
of SM particles in their decays than the squark 2-body decays to neutralinos 
or charginos, particularly if Q — > qxi is the dominant squark decay. On the 
other hand, the LSP in gluino decays must share energy with more particles, 
producing less on the average (see Fig. ||). 

Gluinos and squarks may also be produced at the Tevatron in association 
with charginos or neutralinos (analogous to W and Z + jet production) . These 
processes can be more important than pair production of gluinos and squarks 
if the latter are kinematically limited. Event signatures are similar to Q and 
g production, but possibly with fewer jets, though the events may still pass 
the selection criteria for the squark and gluino searches. For example, Qxf 
production will have one less jet than QQ production, assuming the decay 
Q q'xf or Q ^ gxi 

Promising signatures for squark and gluino production are (i) multiple jets 
and $!rpl3 and (ii) isolated leptons and jets and $!rpEB 

Jets + fJj, 

Both CDF and D0 have performed searches for events with jets and l^rp. This 
signature has significant physics and instrumental backgrounds. The three 
dominant physics backgrounds are (i) Z vv plus jets, (ii) W riy plus 
jets, where the r decays hadronically, and [Hi) ti~* t plus jets, where the r 
decays hadronically. The $rp in leptonic W decays peaks at Mw/'2. — 40 GeV, 
with a long tail at high Iflj. due to off-shell or high-py W^s and energy mismea- 
surements, so a large Ifij- cut is needed to remove these events. Instrumental 
backgrounds come from mismeasured vector boson, top, and QCD multijet 
events. Backgrounds from vector boson production occur for W ev, plus 
jets events when the lepton is lost in a crack or is misidentified as a jet. The 
same problem can occur when the W is produced in a ti event. QCD multijet 
production is a background whfin jet energy mismeasurements cause false Iflj,. 

The D0 Run la analysisLlI searches for events with 3 or more jets and 
Iprp and with 4 or more jets and $rp. The analysis is described in Table ^. The 
resulting mass limits on squarks and gluinos are shown in Fig. ^ (right, the 
plot containing the CDF results also shows the D0 la results) and were set 
using a RIPS model with the following parameters: Mfj±= 500 GeV, tan/3= 
2, fj, = —250 GeV, and = m~. The efficiency and theoretical cross sections 

were calculated using ISAJEli^ assuming 5 flavors of mass degenerate squarks 
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without top squark production and a detector simulation. 

D0 also has a 3-jet analysisEJ based on 79.2 pb^^ of Run lb data. The 
basic requirements are three jets with Et > 25 GeV and a central leading 
jet {\ri\ < 1.1). The $Jrp may be significantly overestimated if the wrong in- 
teraction vertex is used;^ to reduce this effect, the tracks in the leading jet 
are required to point back to the primary vertex. The is required to be 
uncorrelated in (f) with any jet. A cut on the scalar sum of the Et of the 
non-leading jets, called Ht, effectively reduces events from vector boson back- 
grounds. The leading jet is also required to have Et > 115 GeV because the 
only available unbiased sample to study the QCD multijet background had 
this requirement. These cuts are suminarized in Table ^. Vector boson back= 
grounds are estimated using VECBOSjlHl while the tt background uses HERWIGO 
normalized to the EK2 measured tt cross section. The detector simulation is 
based on the GEANlO program. Two techniques were used to calculate the 
QCD multijet background. One compares the opening angle between the two 
leading jets and the the signal sample to the distribution in a generic 

multi-jet sample. The other selects events from a single jet trigger which pass 
all the selection criteria except for the $!t requirement. The ]^t distribution is 
fit in the low ^j, region, and extrapolated into the signal region. The complete 
set of background estimates can be found in Table 

The D0 data have been analyzed in the context of a minimal SUGRA 
model. For fixed tan/3, Aq, and sign of fi, exclusion curves are plotted in 
the mo — mi/2 plane, Fig. ^ (left). The limits are from_Lhe 3-jet, 79.2pb~"'^, 
analysis only. Efficiencies are calculated using ISAJElEil for production of 
gluinos and five flavors of squarks without stop squark production. For each 
point in the limit plane, the $It ^^'^ cuts are reoptimized based on the 
predicted background and SUSY signal. Figure ^ shows the ]^t ^ function of 
mo and mi/2 for tan/3=2, Aq=0, /i < 0. When mo3>m2/2: the $!t signature is 
degraded, because m~ ^ Mg and thus higher multiplicity gg events dominate. 
Since higher multiplicity also means higher Ht, varying the cuts cap-maintain 
sensitivity. These results are robust within the SUGRA frameworkc^ 

The CDF analysis of the Run lb data set is nat yet complete, but the 
Run la result based on 19 pb~^ has been pubhshedES The basic requirements 
are 3 or 4 jets and 60 GeV of ]^t- The full set of cuts is fisted in Table |[ 
As in the D0 search, the direction of the $!t is not allowed to coincide with 
that of a jet, and events with leptons are rejected to reduce the background 
from W and top events. The variable S, which indicates the significance of the 
$Ti is used to reduce fake ]^t measurements; S is calculated by dividing the 
$Jt by the square root of the scalar sum of the Et in the calorimeters. The 

"^The calculation of I^t uses the event vertex to calculate Et for all objects. 
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Table 4; Selection criteria for Tevatron squark and gluino searches in the 3 or 4 jets+^j, chan- 



nels. Cuts specific to the D0 4— jet analysis are in parentheses. 



Quantity 


Experiment 




D0 


CDF 


Trigger 


fT> 40 GeV 


$Jt>35 GeV 






1 jet with > 50 GeV 




> 75 - 100, (65) GeV 


> 60 GeV with S >2.2 GeV^/^ 




> 25, (20) GeV 


> 15 GeV, \t]\ < 2.4 


leading Eip 


> 115 GeV |77| < 1.1, {N. A.) 


> 50 GeV 




Ht> 100 - 160, (7V.A) 




A(fii between 


5.7° < A0i < 174.3° 


A(/)j > 30° 


and jet «; 


^(A0i - 180°)2 + A(/)2 > 28.6° 


Acjji < 160° 


i—1 is the leading jet 






Leptons 


Veto all (N.A.) 


Veto all 


Vertices 


confirmed, (Only one) 


Any number 



vector boson backgrounds are estimated using VECBOsllJ norniaJized to the 
CDF Wjj data. Top backgrounds are determined using ISAJEIO normalized 
to the CDF measured top cross section. The QCD background is estimated 
using an independent data sample based on a trigger that required one jet with 
Et > 50 GeV. First all analysis cuts (Table U) are applied to this sample except 
for the S cut, the l^rp cut, and the 3 or 4 jets cut. Next the l^j^ distribution is 
fit and the number of events expected to pass the fj, cut is derived. Finally 
the efficiency of the last three cuts is applied to arrive at the final background 
estimate, shown in Table ^. 

The limits derived from the CDF analysis are shown in Fig.|| (right) within 



the RIPS framework (see Sec. 2.5). In RIPS, a heavy gluino implies a heavy 
X]*, so a light squark {m~ w M~q) decay will not produce much ]^rp. The 

consequence is an apparent hole in the CDF limit for small m~ and large Mg. 
However, lighter gluinos always produce a large because of the enforced 
mass splitting between Mg and M~o . The results of this analysis do not change 

substantially as parameters are varied within the RIPS frameworkll 

The results summarized in Fig. 6 are complemented by the dilepton+^j^ anal- 
ysis shown in Figs. 7 and 9. The limits on the gluino and squark masses in 
each scenario (minimal SUGRA and RIPS) will be discussed below. 
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Table 5: The number of expected and observed events for Tevatron squark and gluino 
searches in the jets+^y channel after performing the cuts in Table W. 





D0 


CDF 


Analysis 


3 jets 


4 jets 


3 or 4 jets 


4 jets 




79.2 


13.5 


19 


19 


Z U, vv 
it 


1.56 ± .67 ± .42 
1.11 ± .83 ± .36 
3.11 ±.17 ±1.35 


4.2 ± 1.2 
1.0 ±0.4 


13.9 ±2.1 ±6.0 
5.0 ±0.9 ±2.7 
4.2 ±0.3 ±0.5 


2.6 ±0.9 ±1.7 
0.4 ±0.2 ±0.4 
2.2 ±0.2 ±0.4 


QCD multijets 


3.54 ±2.64 


1.6 ±0.9 


10.2 ± 10.7 ±4.2 


3.2 ±3.8 ±1.3 


Total Background 


9.3 ±0.8 ±3.3 


6.8 ±2.4 


33.5 ±11 ±16 


8±4±4 


Events Observed 


15 


5 


24 


6 



> 

Siio 

5" 



80 



50 



50 150 300 

Mo (GeV) 

Figure 5: The l^j, c|i8|tribution frmn simulations of squark/gluino events in the D0 detec- 
tor based on ISAJETEJ and GEANTL3 The simulation used SUGRA mass relations assuming 
tan/3=2, Ao=0, and /i < and seven values of mo and mi/2. The numbers in the upper 
right-hand corner of each plot are the mean and RMS of the distribution. The normalization 

is arbitrary. 



Dileptons+^T 

If, in the cascade decay chain of the Q's and g's, two charginos decay xt ~^ 
(vXi: or one neutrajino decays X2 ~^ ^^^~Xij the final state can contain 2 lep- 
tons, jets, and $!rpEB This channel has the advantage of being relatively clean 
experimentally. The requirement of two leptons significantly reduces jet back- 
grounds and removes most of the W backgrounds. Requiring that the mass of 
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CDF Preliminary 




Figure 6: (Left) The D0 Excluded region in the mo — va-^j^ plane with fixed parameters 
tan/3 = 2, Aq = 0, and < 0. The heavy solid line is the limit contour of the D0 jets and 
missing transverse energy analysis. The dashed line is the limit contour of the D0 dielectron 
analysis. The lower hashed area is a region where mSUGRA does not predict EWSB cor- 
rectly. The hashed region above is where the sneutrino is thei^SP. (Right) The CDF mass 
limits on squarks and gluinos from the search in jets and ^yPi using 19 pb~^ of data and 
the ISAJET 7.06 Run I Parameter Set (RIPS) with the indicated values (solid area). For 
77i~ < Mg, the cross section used is leading order, and three or more jets are required. For 

Mg < the cross section is NLO,El and four jets are required. The line labelled "D0 

PRL" is the D0 result from Run la using 13.5 pb-i of datai] 



the two Icptons be inconsistent with the Z mass removes most of the rest of the 
vector boson backgrounds. If the leptons are required to have pt>20 GeV, the 
major background from physics processes is tt^ bW^bW~ bb£^£~]^rp. As 
the cut on lepton pT is lowered, Z — s- r+r", where the r's decay semileptoni- 
cally, also becomes an important background. The instrumental backgrounds 
are small. The spectacular signature of like-sign, isolated dileptons, which is 
difficult to produce in the SM, can occur whenever a gluino is produced directly 
or in a cascade decay, since the gluino is a Majorana particle. This property 
is exploited in the CDF dilepton searches. r— in 

Figures ^ and ^ show the 000 and CDFE^'O results from Run lb, once 
again compared to SIJGRA and RIPS, respectively. The CDF limit is based on 
NLO cross sectionsp^ and the D0 limit on LO cross sections. The D0 limits on 
niQ and mi/2 are calculated including contributions from the production of all 
sparticles (for instance, associated production of neutralinos or charginos with 
squarks or gluinos), while the CDF result only considers Q and g production. 
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Table 6: Selection criteria for Tevatron searches for squarks or gluinos in the dilcptons, 2 

jets and channel. 



Quantity 


Experiment 


D0 


CDF 


jCdt (pb-^) 


92.9 


81 


J-Jrp J ±_Jrp 


> 15, 15 GcV 


> 11,5 GcV 


Pt ^Pt 


N.A. 


> 11,5 GcV 




> 20 GcV 


> 15 GeV 


Mass window cut 


Mz ± 12 GcV 


N.A. 


AR between leptons and jets 


N.A. 


> .7 




> 25 GcV 


> 25 GeV 


like sign dilcptons 


no 


yes 



Tabic |6| gives the selection criteria for the two analyses. The experimental cuts 
axe chosen to identify two high-pT leptons, which come predominantly from 
Q and g decays into charginos or neutralinos which in turn decay into real or 
virtual W m Z bosons. 




Figure 7: (Left) The p0 limits on the SUGRA parameters mo and miy2 from the 2 leptons, 
2 jets, and searchEj for tan/3=2, Ao=0, and < 0. (Right) The same plot for tan/3=6, 
Ao=0, and fi < 0. In both plots, the dark shaded area is the region in which SUGRA does 
not produce electroweak symmetry breaking. Selected contours of squark and gluino mass 

are also shown. 
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CDF Preliminary 




50 100 150 200 250 300 350 400 450 500 

gluino nass (GcV/c^) 



Figure 8: CDF limits on the squark and gluino masses from the 2 like sign leptons, 2 jets, 
and If^rj. search in 81 pb~^. The limits were set using the ISAJET 7.06 Run I Parameter Set 
(RIPS) with the indicated values. 



D0 has also presented an experimental limit in the Mg — m~ plane Fig. ^ 
which allows a comparison with the CDF limit Fig. ||. For m~ ^ Mg or, 
equivalently, for mo 3> toi/2, 99 pair production is the dominant SUSY process. 
As mo(TO~) is varied with the other parameters fixed, the branching ratios for 
the 3-body gluino decays to charginos or neutralinos and jets become fairly 
constant, so the production rate of leptonic final states becomes constant; the 
experimental limit approaches a constant value asymptotically, as can be seen 
in both the D0 and CDF plots shown in Figs. || and |9[ respectively. Observe 
that, for large enough values of the gluino mass, the leptons easily pass the 
experimental cuts, so the experimental efhciency also becomes constant. 

The relation m~ ^ Mg is not possible in SUGRA, and is treated in an 
ad hoc manner in RIPS. There is no limit in this region for either opposite- 
or like-sign dilepton pairs because the large, fixed slepton masses limit the 
branching ratios to leptonic final states. The possibility of like-sign dilepton 
pairs is further reduced because both the gg and gQ cross sections (which 
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produce like-sign leptons because the gluino is a Majorana particle) and the 

QQ cross section (which produces like-sign leptons because the squarks have 

the same charge) are small in this region. It is very difficult for QQ* production 

to yield like-sign leptons in general. 

When m~ ~ Mg, the gg cross section is supplemented by the gQ cross 

section. Just above the diagonal line at Mg = m~ {i.e. m~ just larger than 

n n Q Q 

Mg) in Figs, g and there are "noses" in the limit plots, with the limit 

becoming stronger close to the diagonal. 




50 100 150 200 250 300 350 400 
Gluino mass (GeV/c^) 



Figure 9: Excluded region from various D0 analyses in the m Mg plane with fixed 

Q 

mSUGRA parameters tan /3 = 2, j4o=0, and < 0. Note that there are no mSUGRA 
models in the region to the right of the diagonal thin line. The heavy solid line is the limit 
contour of the D0 Run lb 3 jets and missing transverse energy analysis. The dashed line 
is the limit contour of the D0 Run lb dielectron analysis. The dot-dashed line is the limit 
contour of the D0 Run la 3 and 4 jets and missing transverse energy analysis shown only 
in the region with valid mSUGRA models. 

The limits in Figs. ^ and ^ are for a specific choice of parameters within the 
RIPS or SUGRA framework. If /z, At and tan/3 are varied, the branching ratios 
into charginos or neutralinos can vary strongly. The sensitive dependence on 
the parameters can be seen within minimal SUGRA models from the D0 limits 
in Fig. 0. The dip in the tan/3=2 limit (left), around niQ— 70 GeV, is a 
point where > M~o > ma and BR(x2 ~^ ^'^Xi) — 1, so the detection 
efficiency is very sensitive to the choice of high energy parameters mo and 
mi/2- In Fig- (right), with tan/3=6, the limits are severely reduced compared 
to Fig. 1^ (left), with tan/3=2, in the region where the squark mass is large 
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compared to the gluino mass. For large tan /3, the mass splitting M~± ~o — M~o 

is reduced, so that the leptons from the Xi and X2 decays are softer. The 
non-trivial shape of the limit curves results from an interplay between the 
cross section being larger when mo and mi/2 are smaller (sparticle masses are 
smaller) and the mass splittings being smaller. Consequently, although the 
dileptons+jets+^-p signature is an excellent discovery channel with little SM 
background, it is hard to set significant parameter limits even using SUGRA 
models. 

From the present analyses in the ^j^+jets and dileptons+^^j- channels, 
some preliminary conclusions can be drawn on the squark and gluino masses. 
These depend, however, on the assumed SUSY parameters. The D0 limit on 
the gluino mass effectively develops a plateau for large mo at 185 GeV for 
tan/3= 2, and at 134 GeV for tan/3= 6. The CDF limit on the gluino mass 
is 180 GeV for tan f3= 4 for large m~. Instead, for equal squark and gluino 
masses, the D0 mass limit for tan /3=2 is 267 GeV, using all SUSY production 
and decay modes in the model. From the CDF analyses and m~ ~ Mg, the 
limit is about 220 GeV for tan /3 = 4. A direct comparison of all the above 
results is rather difhcult since D0 and CDF have done analyses assuming 
different sets of MSSM parameters (see Figs. 6-9). Moreover, CDF considers 
only squark and gluino production, while D0 considers all possible sparticle 
production, and the associated production of neutralinos or charginos with 
squarks or gluinos can have an impact on the experimental limits. 

It would be very useful for purposes of comparing and combining the two 
experimental limits to have both collaborations use at least one common model 
(such as SUGRA), and agree on several values of the parameters to do the 
searches. For example, the two collaborations could present their limits in 
the mi/2 — tan/3 plane (for large mo). Secondly, they could move (partially) 
towards more experimentally-based quantities by plotting contours of cross 
section limit and also contours of acceptance x efficiency in the mo-mi/2 plane. 
This would eliminate the strong model dependence on the branching ratios. 
The experimental acceptance for the signature of two leptons+jets is much less 
model-dependent, since it simply reflects the hard kinematics from the decays 
of two heavy objects. A presentation of cross section x branching ratio limits, 
in addition to the mass limits, would be of more general use to model-builders. 

4-3 Top Squarks 

The top squark (stop) is a special case worth a separate discussion00 The 
mass degeneracy in the stop sector is expected to be strongly broken, and, 
for sufficiently large mixing, the lightest stop can be expected to be rather 
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light, possibly lighter than the lightest chargino. The lightest stop has about 
a tenth the production cross section of a top quark of the same mass, because 
the threshold behavior is (compared to (3 for fermion pairs) and only half 
the scalar partners are being considered. At leading order, the cross section is 
independent of the gluino mass and depends only on the stop mass^Q Due to 
the large left-right mixing, the NLO SUSY QCD corrections must deal with 
different left- and right-handed couplings of the quarks to squarks and gluinos. 
The results for the stop-pairjDroduction cross section as a function of the stop 
mass are plotted in Fig. [l^.L3 



stop and Top Pair Production 




80 100 120 140 160 180 200 



(GeV) 

Figure 10: The cross section for pair production of one stop quark versus the stop mass at 
the Tevatron, calculated to NLO. The band represents the variation of the scale from 1/2 to 
2 times the stop mass. The top pair production cross section versus the top mass is shown 
for reference, calculated only at LO. The band represents the variation of the scale from 1/2 
to 2 times the transverse mass pT © mt , where px is the transverse momentum of the top 

quark. 

The stop caxL be produced directly as tt* pairs or, depending on the stop 
mass, indirectlyllj in decays of the top t ix\^ or sparticles, such as ~* bt. 
Also depending on the stop mass, one of three decay modes is expected to 
dominate. If (a) m^^ > m~± + m^, then ii can decay into hxi , followed by 
the decay of the chargino. This can look similar to the top decay t bW , but 

^ Since the incoming partons are not tops, then if flavor changing neutral currents are 
suppressed by a Supersymmetric GIM mechanism, stop production via gluino interchange is 
not allowed. After including the NLO SUSY QCD corrections the dependence on rrig and 
stop mixing becomes explicit but, in practice, numerical results are insensitive to the exact 
gluino mass and the mixing. 
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with different kinematics and branching ratios for the final state, fnstcad, if the 
stop is the Hghtest charged SUSY particle, it is expected to decay exclusively 
through a chargino-bottom loop as (b) ti —f cx^, which looks quite different 
from SM top decays. Finally, the stop can decay (c) i — > bWxi or if the stop 
is quite heavy into t txi 

The possible signals from it* production, with decay (o) and depending on 
the chargino decay modes, are: {i) bbl^i~$rp^ (ii) bW^jj$rp^ or {in) bbjjjj$rp. 
These are similar to ti final states, except (iii) has real If decay (6) 

dominates, this yields a signature of 2 acoUinear charm jets and $]rp. Finally, if 
decay (c) occurs, the events are similar to ti events, except that the kinematics 
of the individual t and t are altered and there can be much more l^j.. 

If the stop is in the range 100 — 150 GeV, ii* production may be too small 
to observe and it might be easier to observe a light stop in top quark decays 
t — > ixi- If Xi is Higgsino-like, then the BR{t ixi) can be 50%. If decay 
(a) occurs, then top quark events have the same signatures as in the SM but 
they have more $!rp and softer jets and leptons. In case (6), fewer leptons and 
jets are produced and the distribution is affected. If there is one SM top 
decay and one SUSY top decay, the final state can be b£^c]^j,^ which would 
appear at a small rate in the Wjj sample, but not in the SM ti event sample. 

An indirect limit can also be set on the decay (6). Such decays would not 
fall in the SM ti dilepton or lepton+jets samples, but instead would deplete 
them. Given a theoretical prediction for the ti production cross section, the 
branching ratio for decays which deplete the SM ti samples can be bounded in 
a straightforward manner. If decay (a) occurs, the analysis is more involved, 
since the kinematic acceptance for the stop decays must be calculated for many 
different choices of MSSM parameters. Also in case (a), some tt* events will 
feed into the top quark event samples. 

Direct Top Squark Pair Production 

D0 haS|_aearched for ii* productior0 with i cxi using 7.4 pb^^ of Run 
la data.Lj The signature is two acoUinear jets and ^j,, satisfying the selection 
criteria in Table 0. The dijet cross section at the Tevatron is large, and thus this 
signature has large instrumental backgrounds. It also has backgrounds from 
vector boson production. The QCD and vector boson backgrounds would, 
naively, be a factor of l/a^ larger than for the hadronic squark/gluino search, 
as this search requires only 2 jets while the latter searches require at least 3 or 
4 jets. This is not the case since the multijet backgrounds can be controlled 
by requiring A0 > 45° between the l^j. and each jet, and that the jets not be 
back-to-back. The vector boson backgrounds are controlled by requiring that 
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Table 7: Selection criteria for the D0 2 jet+_^j, hadronic direct stop pair production search. 
ji and j2 are the leading and sub-leading jets ranked by Et- 



Quantity 


D0 




> 30 GcV 




90° < A</-< 165° 


A(j> jl-$r^ 


10° < A</« < 125° 




10° < A(/> 


Lepton veto 


e and ^ 


$T 


> 40 GcV 



the two leading jets are separated by at least A(j> > 90°. After these cuts, the 
dominant backgrounds are from W and Z boson production and decay, with 
the largest being W tv. If the r decays hadronically, only 1 additional jet 
is necessary to fake the signature. Top quark production is not as important 
a background for the light stop search as for the conventional hadronic squark 
search because of the lower jet multiplicity requirement. As with the hadronic 
squark/gluino searches, the cuts are not very efficient for signal events. The 
efficiency is largest when the stop is heavy compared to the Xi (near the 
kinematic boundary for the decay ti — > tx\), reaching a maximum value of 
only 4%. The mass difference mi — M~q determines the Et of the charm jet 
and rapidly limits this search mode as the charm jets become too soft (see Fig. 
0). 

With the assumption that BR(ti cx5)=l, the predicted SUSY final 
state depends only on M~o and m^^ . The result of this search is a 95% C.L. 



exclusion limit on a region in the M—q — to^^ plane, shown in Fig. 11, The 



production rate has been calculated using only LO production cross sections 
evaluated at the scale = 2stu/{s'^ + + u^) from ISAJEipil so the limit will 
change somewhat if re-evaluated with NLO production cross sections. 

CDF and EUZl have also presented results from a search for iiil production, 
with ii bxi'B The CDF search is in the lepton+jets chapiiel, and uses a 
shape analysis of the transverse mass|] of the lepton and ^^.£3 The selection 
criteria are given in Table |^. For both searches, the detection efficiency is 
smaller than for tt production with a top quark of the same mass because of 
the softer leptons from the 3-body decay xt ~^ ^^Xi- The mass splitting 



s Transverse mass squared is Af|, = + — {p^ + -?^t)^i ^'^'^ ^ useful ex- 

perimental quantity when information about the longitudinal component of momentum is 
missing. 
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95% Exclusion Limit on t^— >cx° 




50 60 70 80 90 100 110 120 



t, (GeV) 



Figure llr-iMass limits from the D0 search for it* production with the decay i —> at the 
TevatronLJ The decay is kinematically forbidden in the two solid grey regions. The hashed 
regions marked 0j show the LEP excluded regions as a function of the stop mixing angle, 
which determines the strength of the stop coupling to the Z boson. The mixing does not 
affect the tree level process at hadron colliders. 



Table 8: Selection criteria for the CDF lepton+jets+6-tag direct stop pair production search. 

e or fj,. 



Quantity 


CDF cuts 




> 20 GcV 


$T 


> 20 GeV 


E-ij^ , rj range 


> 15 GeV, I77I < 2.0 


E^rji , 77 range 


> 8 GeV, It^I < 2.4 


Jets 


<4, > 8 GeV 


SVX 6-tag 


1, >8 GeV 



Mi 



M~± sets the efficiency for detecting the jets. The results of the CDF 
search are shown in Fig. 12 (left). The decay Xi 



W*Xi is assumed using 



the masses ii) M~± 



80 GeV and M~ 



30 GeV and {ii) M~± 

Xi 



70 GeV 
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Table 9: Selection criteria for the D0 dieloctron+jcts+_^j,direct stop pair production search. 



Quantity 


D0 cuts 


pei Tpe2 


> 16,8 GcV 




> 30 GcV 




< 90 GcV 




< 60 GcV 


tr 


> 22 GcV 




Figure 12: (Left) The CDF cross section limit on direct production of the top squark using 90 
pb~^ of data. The decay mode is f — > 6x]*^(— > VF*Xi)- One W must decay semi-leptonically 
giving a signature of a lepton, ]pj,, and jets. The theoretical cross section is from IS A JET 
7.06. (Right) The D0 95% confidence level cross section limit on the cross section for stop 
production times the branching ratio to p-pnal state containing 2 electrons as a function of 
the mass of the i is shown as a solid lineEj The mass of the lightest charginroi is assumed to 
be 47 GeV. The predicted cross section times branching ratio from ISAJETtll is also shown 

as a dashed line. 



and M~o ~ 30 GeV. Given these mass choices, there is Uttle other parameter 
dependence. Presently, the cross section hmits are above the predicted cross 
sections due to the high Et cuts. 

D0 searches in the dilepton channeH using the cuts hsted in Table ^ 
The signature is similar to the squark and top dilepton searches. The results 
are shown in Fig. Ill (right), assuming M~± = 47 GeV and M-, = 28.5 GeV. 

A substantial background comes from Z — s- r+r^, again requiring a high 
threshold for the Et cuts, and no limit can be set. 

The above analyses were done in regions of SUSY parameters that have 
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been excluded by LEP. They show, however, the procedures to be followed in 
redoing these studies for other regions of the MSSM parameter space. 

Top Squark Production From Top Decays 

CDF has presented another analysis using the SVXptagged lepton+jets sample 
to search for the decay t hXii with ii ^'xf E2l If one of the top quarks 
in a tt event decays t bW(~^ Iv) and the other t iiXi followed by fi 
bxf{-^ jjXi) ort^ 6W(-> jj) and t iiXi followed by ii &xf ^'^Xi), 
the signature is bWvjj+$rp, the same as in the SM, but where the includes 
the momentum of the x?- The lepton+jets channel has a large number of 
events, so a kinematic analysis can be performed on the event sample. Due to 
the mass of the x^ and the intermediate sparticles in the decay chain, the jets 
from the SUSY decay are significantly softer. This difference is exploited as 
the basis of the search. 

Table 10: Selection criteria for the CDF search for top decaying into stop in the signature 
of 1 lepton {£), ^j,, and 3 jets including at least one b tag, where £ = e or fi. The quantity 
I cos 9* I is the polar angle of a jet in the rest frame of the £, fj, and jets. ARi is the distance 
between a jet i and the next nearest jet in r] — <j) space. The jets are ordered in Et, so 

E^> E^> E^. 



Quantity 


CDF cuts 


JCdt 


110 pb-^ 




> 20 GcV 


$T 


> 45 GeV 


Mt{1$t) 


> 40 GeV 


Vt{1$t) 


> 50 GeV 


ll/rp 


> 20 GeV, \ri\ < 2.0 




> 15 GeV, \ri\ < 2.0 


1 cos 6**1 1,2,3 


< 0.9,0.8,0.7 




> 0.9 


Number of SVX 6-tags 


> 1 for E^. >15 GeV 



The cuts listed in Table |lO| are optimized for acceptance of the SUSY decay 
and rejection of M^+jets background. A likelihood function is computed for 
each event reflecting the probability that the jets with the 2"^^ and 3"^ highest 
Et in the event are consistent with the stiffer SM distribution (as compared 
to the SUSY distribution). The distribution of this likelihood function shows 
a significant separation of these two hypotheses. After applying the cuts listed 
in Table O, nine events remain, all of which fall outside of the SUSY signal 
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region. For stop masses between 80 and 150 GeV and chargino masses between 
50 and 135 GeV, a BR(t ^ iiXi)=50% is excluded at the 95% C.L., provided 
that M~o = 20 GeV. Because M~n is fixed in this manner, it is not related to 

M~± as in SUGRA. At present, only this one example is available (for M~o 

already excluded by LEP); more statistics will significantly improve it. 



4-4 Sleptons 

At hadron colliders, sleptons can only be pair produced through their elec- 
troweak couplings to the 7, Z and W bosons. Figure shows the cross sections 
as a function of the corresponding slepton mass compared to the differential 
Drell-Yan pair production cross section, dau^y /dQ, where Q — 2Mx- The 
rate for slepton pair production is at most a few tens or hundreds of fb at the 
Tevatron, and so far neither collaboration has presented results on searches 
for sleptons in the SUGRA or RIPS framework (we describe limits in gauge 
mediated models later). 



Slepton Pair Production at LO 




Mx (GeV) 



Figure 13: Production cross sections for sleptons versus slepton mass at the Tevatron. The 
Drell-Yan cross section for producing e+e" (curve b) is plotted as a function of M^+^- = 
2Mx- On this scale, the Drell— Yan peak (Z — > e+e~) would appear at Mx = 45 GeV. 
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A (stable) charged slepton is not a viable LSP candidate, so the decays 
If ^ ^^Xi* or — s- vxf are expected. The sneutrino, instead, can be the 
LSP, or it can decay invisibly v — *■ vxi^ or visibly i> — -> xf^^- If ™!5 < "z^^ < 
M~o, then the decay £ t'v'v (or ^ qqv) is possible. Promising signatures 

are (i) e^e^, r+r^ plus .^^,0 (m) e^,eT,^T plus .^y, and (iii) e,/i or 
T+jets plus $!rp (or jets plus Although charged slepton production can 

lead to charged leptons in the final state, there is no guarantee. 

The major background to same-flavor lepton pairs is Drell-Yan pair pro- 
duction, with fake from mismeasurement of the lepton or jets in the event. 
Most of this background can be removed by vetoing on a dilepton mass window 
around the Z mass, by requiring significant ^x(> 25 GeV), and by vetoing 
events with the fj- pointing in (j) along one of the leptons or a jet. Top quark 
production is also a major background to a slepton heavier than the lighter 
gaugino, as it produces dilepton events that have real ]^rpE3 Untangling a few 
heavy slepton events from top events would be difficult at the present low level 
of statistics. 

Inclusive searches have the advantage of a larger acceptance than searches 
in exclusive channels. A unique signature of slepton production directly or in 
cascade decays would be the apparent violation of lepton universality. If the 
sleptons are not degenerate, both the production and decays of the sleptons will 
favor one or two leptons over the others, resulting in an imbalance in the de- 
tected e/n/r ratios in SUSY-enhanced channels. The dominant backgrounds 
to inclusive leptons come from hea|5iy flavor production {e.g. 6-quarks), and 
(single) W and Z boson production.Ej Because sparticles are produced in pairs, 
it may be possible to discriminate against SM backgrounds by requiring the 
identification of a part of the decay of the second sparticle. Examples of chan- 
nels that may have enhanced SUSY contributions over SM backgrounds (and 
hence possibly apparent lepton universality violation) are those that have, in 
addition to the lepton, a 7, W, Z, additional lepton, or a third-generation 
particle. 

4.5 Charged Higgs Bosons 

Even though Higgs bosons are not sparticles, the discovery of one or more 
would be considered indirect evidence for SUSY.^ If it is light enough, the 
charged Higgs boson can be produced in the decay of the top quark 
t bH^ SB The branching fraction for this decay depends on the charged 

''If M^± ~ 300 GeV, then there must exist extra hght— matter fields beyond the SM to 
partially cancel the contribution to BR(fe — > S7). 
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Higgs boson mass and tan/3. When kinematically allowed, this branching ra- 
tio is larger than 50% for tan/3 less than approximately 0.7 or greater than 
approximately 50, and completely dominates for very small or very large val- 
ues of tan/3. However, as discussed in section 2.1, values of tan/3 < 0.6 — 0.7 
or above 60 would be associated with large top or bottom Yukawa couplings, 
which will become infinite at scales not far above the TeV scale. In general, at 
reasonably small values of tan/3, the charged Higgs boson decays i/+ cs; at 
large tan/3 it instead decays t^Vt- 

GDF has searched for the decay t — > hH^ using both directH'EzI and in- 
directEil methods. Direct searches look for an excess over SM expectations 
of events with r leptons from the charged Higgs boson decay H'^ t^v^ 
(dominant for large tan/3). On the other hand, indirect searches are "disap- 
pearance" experiments, relying on the fact that decays into the charged Higgs 
boson mode will deplete the SM decays t bW, decreasing the number of 
events in the dilepton and lepton-|-jets channels. 

The CDF direct search at large tan /3 uses two sets of cuts, listed in Table 
pr] , to search for an excess of t's in tt events. The first set selects a sample 
containing a r that decays hadronically, an SVX 6-tagged jet, l^j. and objects 
indicating activity from the second top decay: a second jet and a third jet 
or lepton. As Mh± approaches nit, the b produced in the top decay t hH 
becomes less energetic, causing a reduced efficiency for the jet and 6-tagging 
requirements. To maintain efficiency in this region, a second set of cuts accepts 
events that have two high-i?T t's and 

The signature for hadronically decaying r's is a narrow jet associated with 
one or three tracks with no other tracks nearby. Typically the r tracks are 
required to be within a cone of 10° with no other tracks within a cone of 30°. 
Fake rates, measured as the probability that a generic jet is identified as a r, 
are approximately 1% or less. These fake rates are too high to identify r's in 
a sample dominated by QCD. However, if another selection criterion is added 
that further purifies the sample, r's can be identified with a good signal-to- 
background ratio. For example, hadronic decays of the r are observed in (i) 
monojet {W Ti/)jr-lii^ lepton, and jet {Z tt) and {in) the Icpton+jets 
top quark samples .E§E2l 

The CDF direct search for tt events with one or two charged Higgs bosons 
decaying to r leptons sets limits by two methods. In the first method, a 
tt production cross section and a SUSY model {Mjj± and tan (3) is assumed 
and the number of expected r events is computed. If the number expected 
is too large to be consistent with the observed number at the 95% C.L., the 
SUSY model is excluded for that ti cross section. This method excludes a 
charged Higgs with mass less than 155 GeV (100 GeV) if tan/3 is greater than 
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approximately 100 (50) and the top cross sectionj is 7.5 pb, as shown in Fig. 14 
(left). 

The second method combines the observation of i ^ bW decays into lep- 
tons (e or fi) and jets with the number of r decays from the direct search. 
This has the advantage that a top production cross section does not need to 
be assumed. The lepton+jets sample defines a top production cross section 
which, in turn, through the SUSY model, predicts the number of r events 
expected. If the number is too large to be consistent with the observation, 
the modeLis excluded. The limits set by this method are presented in detail 
elsewhere.Ell Qualitatively, the limits are similar to those set by the indirect 
method (discussed below). 



Table 11: Selection criteria for the CDF direct search for t — > bH^{-^ tu) in 100 pb~^ of 
data. The r's are identified in their hadronic decay modes as one or three isolated, high-py 
tracks. Events are accepted if they pass the cuts in either analysis path. 



Quantity 


CDF 


Analysis path 1: 




> 20 GeV 


tr 


> 30 GeV 


£;fi , SVX tagged 


> 15 GeV 




> 10 GeV 


Additional object 


e, /i, T or 




3''^ jet with Et > 10 GeV 


Analysis path 2: 




> 30 GeV 


A0(ri,T2) 


< 160° 




> 30 GcV 



At small tan/3 a direct search is difficult since the charged Higgs decays 
into two jets. Instead only the indirect method is applied. 

The indirect method can be applied to both small and large tan/3 searches. 
The observed numbers of dilepton and lepton+jets events are consistent (at 

* The reader should be aware that there are subtleties in analyses that assume cross sec- 
tions. The CDF experiment normalizes all cross sections to its measured proton-antiproton 
cross sections, rather than to a hard process such as W production. The CDF total cross 
section, in the judgement of one of the authors (HF), is most likely 10% too high, and thus all 
cross sections (in particular the top cross section) are too high. For analysis which compare 
different channels internally this has no effect, but for analyses which compare to theoretical 
predictions the reader should be careful. The DO experiment normalizes to a weighted mean 
of the CDF and E710 values for the total cross section which is 2.4% lower than the CDF 
value. 
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the 95% C.L.) with a minimum production and decay rate in the SM. For 
an assumed tt cross section (and a charged Higgs boson mass and tan/3), a 
simulation of the expected mixture of SM and charged Higgs boson decays 
predicts a number of events in the SM channels. At points in the parameter 
space where the decays to charged Higgs bosons are more pronounced, the SM 
contributions are diminished. If the number expected in a SM channel is not 
consistent with the rate defined by the data, the assumed values are excluded. 
This method provides the limit displayed in Fig. ^ (left). The cross section of 
5 pb is the expected cross section for a top mass of 175 GeV; the curves using 
7.5 pb show the sensitivity of the limit to the assumed top cross section. Also 
shown in the figure is how the limit in the region of large tan f3 can be extended 
using the assumptions of the indirect method. In this region the possibility 
that a r decay produces a high-p^ lepton is included. 

The area in Fig. |4| (left) labeled "ratio method" is the exclusion region for 
an indirect search that does not make an assumption for the ti cross section. 
If charged Higgs boson decays were competing with SM decays, the ratio of 
dilepton events to lepton+jets events would decrease, regardless of the ti cross 
section. This occurs because the lepton+jets yield is proportional to the SM 
branching ratio while the dilepton yield is proportional to the SM branching 
ratio squared. For each SUSY parameter point, the lepton+jets sample can be 
used to infer a top cross section which, in turn, predicts a number of dilepton 
events. The point is excluded if the prediction is inconsistent with the dilepton 
data. Although this method excludes less parameter space, itis,important since 
the tt cross section may be enhanced by SUSY mechanismslllS such as g ^ tt. 
At present, this method only excludes values of tan/3 ~ 0.7, which are not of 
much interest according to present theoretical bias. 

D0 has also searched for a charged Higgs boson lighter than the top quark 
using the indirect niethodE3 The analysis compares the number of events ob- 
served in the lepton+jets channel to the number predicted assuming a theoret- 
ical tt production cross section. The limits depend on the mass of the charged 
Higgs, tan/3, and the top quark mass mt. Table shows the selection criteria 
used in the search. Fig. |lj (right) shows the excluded region. 

Recent studies have shown that quantum SUSY effects (SUSY QCD and 
electroweak radiative corrections) to the decay mode t bH~^ (with subse- 
quent decays into r's) may be important and should be considered in future 
analyses .Ell 
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Figure 14: (Left) Exclusion space for the CDF searches for charged Higgs boson decays of 
the top quark in ti events. The shaded regions are from the indirect searches. For the 
regions labeled o-jj = 5.0 and 7.5 pb a top production cross section is assumed and points 
are excluded if the predicted SUSY decays have depleted the SM channels to an extent that 
they are inconsistent with the data. The "ratio method" is an indirect method comparing 
the number of lepton+jets events to the number of dilepton events and no top cross section 
is assumed. The region excluded with solid lines at high tan /3 is from a direct search for 
events where one or both top quarks in a tt event decay to bH^{-^ t+i/) and information 
from the SM channels is ignored. (Right) The results of a D0 indirect search for a charged 
Higgs boson assuming mt = 175 GeV and a tt production cross section of 5.53 pb and 4.77 
pb. This limit is based on the full Run I D0 data sample. 



^.6 Neutral Higgs Bosons 

Within the MSSM, the main production channels for the hghtest CP-even 
Higgs boson h at the Tevatron are the same as for a SM Higgs boson, Wh or 
Zh production.^ The cross sections behave in such a way that these channels 
are relevant for large values of the CP-odd mass Ma (the SM limit) or for 
small Ma and small tan j3. The heavy CP-even Higgs boson H could become 
marginally relevant for searches at an upgraded Tevatron through ZH , WH 
production, in some restricted region of parameter space, complementary to 
the one relevant for the light CP-even Higgs boson searches. In addition, the 
enhancement of the bottom Yukawa coupling in the large tan [3 regime can 
render the production processes hbb—Ahh, and Hhb useful to perform searches 
in a large region of parameter space J23 

Both collaborations have searched for a neutral Higgs boson in the mode 
qq' ^ W* ^ W{-^ ev,iiu)h{-^ bb). D0 has searched in 100 pb'^ of data 
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Table 12: Selection criteria of the D0 search for a charged Higgs boson produced in top 
quark decays. In addition, events are vetoed if the is aligned in ip within 25° of a muon, 
or if the muons in a event with a ^-tagged jet have a good fit to the decay Z — > fifi. 



Quantity 


D0 topological 


D0 tagged 


Et threshold on leptons 


20 GeV 


20 GeV 


Max 77 for leptons 


2 (e) 1.7 


2 (e) 1.7 (m) 


Number of jets 


4 


3 


jet Et threshold 


15 GeV 


20 GeV 




25 GeV (e) 20 (^i) 


20 GeV 


Hx 


180 GeV 


110 GeV 


Sum of lepton Et and ]frp 


60 GeV 


N.A. 


Aplanarity 


0.065 


0.04 


\vw\ 


2.0 


N.A. 


^-tagged jets 


veto 


require 



using a data sample containing a lepton, ]/]rp and two jetsEi One of the jets 
must have a muon associated with it for 6-tagging. The cuts are listed in 
Table 13, Twenty-seven events pass the selection criteria; 25.5 ± 3 events are 
expected from Wjj and ti. The limits shown in Fig. |l^ are set by a simple 
event-counting method and by fitting the bb dijet mass spectrum. 

CDF has receittly completed a similar search for the same decay mode using 
109 pb"-'^ of dataEj All events must have one SVX 6-tag. These events are split 
into single-tagged (one SVX tag) and double-tagged samples (two SVX tags 
or one SVX and one lepton (e or fi) tag). The 36 (6) single-tagged (double- 
tagged) events are consistent with the 30 ± 5 (3.0 ± 0.6) expected from SM 
T4^+jets and tt. Both the single- and double-tagged dijet mass distributions 
are fit simultaneously to set the limits shown in Fig. |l^. 

The process qq ^ Z* ^ Zh occurs at a comparable rate to the W* 
process. CDF has searched for both associated production processes assuming 
W/Z ij Ej The event selection criteria are listed in Table |3[ In 91 pb^^ of 
data, 589 events remain, consistent with the expectation from QCD heavy- 
flavor production and fake tags. To set limits, the bb dijet mass spectrum is fit. 
Also shown in Fig. |l^ is the SM production cross section for Wh and Zh as a 
function of the Higgs boson mass. The present experimental limits are roughly 
two orders of magnitude away from the predicted cross section. However, Run 
II will provide at least 20 times the data. This plus the possibilities of looking 
at other decay modes (i.e. Zh — s- vvbb) holds promise for Higgs physics at the 
TevatronEj 
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Table 13: Selection criteria of Tevatron searches for the associated production of a neutral 



Higgs boson and a W or Z, and the Higgs boson decays to bb. 



Quantity 


CDF 


D0 


WH ^ iiybb 




> 25(20) GeV 


> 20(20) GeV 




> 25(20) GeV 


> 20(20) GeV 




> 15 GeV 


> 15 GeV 


^-tagging 


one SVX tag 


one /i tag 


{W,Z)H^jjbb 


Quantity 


CDF 


rjrp 


> 15 GeV 


6-tagging 


2 SVX tags 


Pribb) 


> 50 GeV 



CDF & DO Preliminary 

<t(pp ^ W/Z + X°) X Br(X° ^ bb) 

95% CL upper limit 




_ CDF (W/Z ^ qq) - 90.6 pb - 
_ CDF(W^lv)-109pb"' 
. . DO shapes (W -> Iv) 
.... DO counting (W -> Iv) 
DO results from 1 00 pb'^ 




90 100 110 120 130 140 

Higgs Mass (GeV/c^) 



Figure 15: Limits from CDF and D0 for the associated production of a neutral Higgs boson 
and a, W or Z boson. The CDF limits are shown for the final states of £ubb and jjbb, and 
the D0 limit is for the final state £i/bb. The limit is set using a simple counting method and 
by fitting the bb spectrum ("shapes"). 



D0 has also searclied for a fermio-phobic Higgs, i.e. one with suppressed 
couphngs to fermionsEZl For a hght neutral Higgs boson, the decay through 
a virtual W loop to a 77 final state can be dominantEj Events are selected 
containing two photons with Et> 20 and 15 GeV, and two jets with Et> 20 
and 15 GeV. No evidence of a resonance is seen in the mass distribution of 
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the 2 photons, and D0 excludes, at a 95% C.L., such a Higgs with masses less 
than 81 GeV. The branching fraction for /i ^ 77 is taken fr'om Ref. |98|. 



Table 14: Selection criteria for the DO search for a Higgs boson produced in association with 



a hadronically-dccaying W and which decays to two photons. 



Quantity 


D0 


Et threshold on photons 


1 above 20 GeV, 1 above 15 


I77I on photons 


<1.1 or 1.5< It^I <2.25 


Et threshold jets 


1 above 20 GeV, 1 above 15 


I77I on jets 


< 2 


vector sum of photon Et 


< 10 GeV 


vector sum of jet Et 


< 10 GeV 



4-7 R-Parity Violation and a Short-Lived LSP 

Allowing for R~parity violation in the MSSM opens a host of possibilities at 
the Tevatron. Both baryon-number-violating operators (UDD) and lepton- 
number-violating ones [LLE and LQD) are possible. There are many resonant 
and non-resonant particle production mechanisrns_and subsequent decay pro- 



cesses which have been analyzed in the literature.EH3 In this section, we restrict 
ourselves to the experimental analyses performed so far. 

The possible excess of HERA events at large has triggered interest in 
studying the consequences of the interaction of a light_aauark (preferably a 



top or charm squark) with an electron and a d quark.t£3 If the gluino were 
heavier than this squark, then gluino pair production at the Tevatron and 
the decay g — s- ccl through R-conserving couplings, followed by the RPV 
decay cl — > e^d, would yield the signature of two electrons and 4 jets. If 
the RPV decay cl e~^d is allowed through the coupling X'121, then from 
the structure of the R-parity violating Lagrangian (Eq. (p4|)) it follows that 
sl Ved, dji — !■ e~c, and dj^ — > VeS are also allowed. If mg^ ~ to^^ (which is 
guaranteed) ~ mj^ (which is probable), then the gluino decays equally to c^c, 

SiS, and dRd (+ h.c.) final states. Assuming that only RPV decays occur, 
then 1/2 of gluino decays produce a charged lepton. Therefore, gg production 
produces like-sign dileptons 1/8 of the time. The requirement of only RPV 
decays demands M-n > m~. 



CDF has performed a searchllHj considering the RPV squark decays with 
the signature of two like-sign electrons and two jets. In 105 pb~^ of Run la 
and lb data, no events remain after all cuts are applied (see Table |5h. Varying 
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Quantity 


CDF 


J-^rp ^ Jjjrp 


> 15 GeV, \t)\ <1.1 




±2 


j'p ^ Jr-p 


> 15 GeV, |r;| <2.4 


S =]pj, / ^Y,Et 


< 5 GeVi/^ 



Table 15: Selection criteria of the CDF search for R-parity violating processes using 

105 pb-i of data. 



the masses of the SUSY particles does not aher the acceptance significantly 
since they are heavy enough for the decay products to easily pass the Et 
thresholds. Because of this, the limit on the cross section times branching 
ratio is approximately constant at 0.19 pb. For mg^ = 200 GeV, this excludes 
Mg < 230 GeV, assuming BR{gg e±e±X) = 1/8. 

Allowing for possible R-parity conserving squark decays, the decay Q 
qxf is possible, where Xi is the LSP. Since the LSP has no R-parity conserving 
decays kinematically accessible, the R-parity violating decay Xi ~^ cde~ or 
cde"*" occurs through a virtual charm or down squark, while x\ ~^ dsv or dsv 
occurs through a virtual strange or down squark. The exact branching ratio for 
^ + X depends on sparticle masses and the mixing of the neutralinos. 
For the analysis, five squark masses are assumed to be degenerate and any 
squark pair can lead to like-sign dielectron events, since Xi is a Majorana 
particle. Squark masses less than 210 GeV are excluded if the mass of the Xi 
is more than half of the squark mass and the gluino is heavy. For lighter x^j 
the 3-body decay of the x? can produce electrons that are too soft to satisfy 
selection criteria. 



4-. 8 R~Parity Violation and Long-Lived Heavy Charged Sparticles 

If R-parity is violated, and the LSP is charged, it can manifest itself as a long- 



lived charged particle (see Sec. 2^) in a collider detector. The particle can be 
identified by measuring the dE/dx energy loss as it passes through the CDF 
SVX and CTC detectors. For a given momentum, a heavy particle has a slower 
velocity and hence a greater energy loss than a relativistic particle (/3 ~ 1). If 
the particle is weakly interacting or massive enough to kinematically suppress 
showering, it will penetrate the detectors and be triggered on and reconstructed 
as a muon with too much_eiiergy loss. A result using part of the Run I data has 



been presented by CDFE£3 and is updated with the full data set here. In 90 
pb~^ of inclusive muon triggers {pt >30 GeV), CDF searches for particles with 
ionization consistent with /?7 < 0.6 and finds 12 events depositing more than 
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twice the energy expected from a minimum ionizing muon. This is consistent 
with the number of events expected from muons which overlap with other 
tracks to fake a large dE/dx signal. 

The CDF result can be used to exclude some SUSY scenarios with R- 
parity violation (RPV). For example, the lightest tau slepton could be the LSP. 
Its production rate through R-parity-conserving couplings can be determined 
from Fig. [l^. If A333 is the only large RPV coupling, the decay f rvr can 
occur with a hfetime fixed by A333 and ruf (see Sec. ^). For small enough 
-'^333 1 this decay can occur outside the tracker, leading to the desired signal if 
the f is travelling slowly enough. 

4.9 Photon and l^j, Signatures 

An Unexpected Turn: the CDF ee^^$j. Event 

Supersymmetry has so many parameters that the full range of its allowed sig- 
natures may be hard to predict. In ApipiJ995, the CDF experiment recorded 
an event with a very unusual topologyli23 which may have SUSY interpreta- 
tions. It has four electromagnetic clusters, which pass the typical cuts for two 
electrons and two photons, and A display of the event is shown in Fig. |l6|. 

The electron in the central region of the detector is well-isolated and is 
associated with a track that has a pT in good agreement with the e~ hypothesis. 
The two photons are also well-isolated and have no associated tracks. The 
"electron" at large 77 is more difficult to identify positively. The associated 
track should only cross a part of the inner CTC where the occupancy is too 
high to find the track. Hence, its charge cannot be determined. The VTX, 
a wire chamber surrounding the SVX but inside the CTC, measuring in the 
r — z view, has a track at the correct rj for the electron hypothesis. The path 
through the cluster and the event vertex can be searched for tracks in the 
SVX, and this analysis is underway. The probability that the event could be 
produced in the SM, including the probability that one or more of the objects 
is fake, is being estimated. The preliminary results indicate that the number of 
expected ££"f"f$!j, events is many orders of magnitude less than one. However, 
the data set was derived from over three trillion collisions, and the probability 
of all signatures which would be considered "rare" must be estimated (an 
impossible task) to determine the significance of one event. 

There have been two main proposals for a possible SUSY explanation of 
the event: the Gravitino LSP and the Higgsino LSP model (for non-SUSY 



explanations, see Refs. [107|, for example). Both proposals also suggest other 
signatures that should be expected within these models and which are pre- 
sented in the following. The Tevatron collaborations have completed some of 
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Run 68739 Event 257646 
28 Apr. 1995, 22:41:20 



Figure 16; The very unusual CDF event containing two 'electrons', two 'photons' and missing 
By. The display is the calorimeter cylinder unrolled into a plane. The towers represent 
energy deposition, with the height of the tower proportional to Et- 



these searches, which are also discussed below. 



Gauge Mediated Low Energy SUSY Breaking: Gravitino LSP 



The CDF analysis of 
SUSY breaking niodelsli°^*°^*^° 



aVinvp nrip event reminded theorists of low-energy 



113 which had long ago lost favor to SUGRA 



models. In these models the (usually ignored) gravitino (G) is very light and 
becomes the LSP. The lightest SM superpartner becomes the next-to-lightest 
supersymmetric particle (NLSP), which is unstable and xbcays into its SM 
partner plus the Goldstino component of the gravitino.Eij In the simplest 
gauge-mediated models, the squarks are heavy and the gauginos obey the 
unification relationship in Eq. ( ITTI ) . Generically the NLSP can be a neutralino 
or a slepton (most plausibly a right-handed slepton and, due to the larger 
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Yukawa coupling, a f). If the scale of SUSY breaking is not far above the 
electroweak scale (< a few 1000 TeV), the NLSP will decay within the de- 
tector, leading to distinctive signatures as displaced vertices or heavy charged 
sipptnns Hpraying into leptons, possibly with a kink to a minimum ionizing 

If a gaugino-like neutralino is the NLSP, the only modification to SUGRA 
phenomenology, where all sparticles decay down to Xii is that Xi then decays to 
a photon and $!rp. The production of X2xf^ xtxi and £^1]^ pairs, followed by 
cascade decays, leads to the final states WZj"f+$]rj^^Wi'^£~j"f+]^j^, WWjj+ 
If^rp and £+£"77 + ^5-, all with comparable ratestiJ A logical starting place 



the CDF event can be interpretfijd as either ee* production followed by 



for searches is in the inclusive two photon and channdliiSjji particular, 

108 113 



eXi or Xi'Xi production^j followed by ^ e i^eXi- If the coupling 



between the gravitino and matter is large enough, then the lightest neutralino 
can decay xH iG inside a collider detector, yielding the desired signature of 
e^e~^^]prp. However, it follows that if one adjusts the parameters of the model 
to explain the multilepton plus photons CDF event, then a very largtj rate of 
multijet plus multileptons plus photon(s) events is to be expected.lliallia The 
fact than none of these other signatures has been detected makes the above LSP 
G explanation of the CDF 6677^7^ event unlikely. Other possible explanations 
may, however, remain open. 



Signatures of photons-|--^rp can point towards models of low-enn 



SUSY 

breaking, but there are other possible signatures in these models. tillEi3 If the 
NLSP is a neutralino which is mainly Higgsino-like, then x\ decays to the light- 
est Higgs boson (or the heavy CP-even or the CP-odd neutral Higgs bosons if 
they are sufficiently light) plus a gravitino. The Higgs boson will subsequently 
decay into bh. Hence, the signature of 4 6-jets, which reconstruct the lightest 
Higgs boson mass in pairs, plus $rp is possible. If the NLSP is a right-handed 
slepton, then the decay I £G occurs, yielding lepton pairs and ]^rp as fi- 
nal signature of slepton pair production. The dilepton signature will suffer 
from large irreducible backgrounds, but the production and decay of heavier 
sparticles can give spectacular signals. For example, the pair production of a 
left-handed slepton which cascade decays into a right-handed slepton and a 
neutralino can yield six leptons+^^j^ in the final state. Also, since the NLSP 
slepton can be tr, signatures with many r leptons are possible. 

If any of the above signatures were observed experimentally, a measure- 
ment of the decay length of the NLSP would provide information about the 
scale of supersymmetry breaking. However, the scale of SUSY breaking might 
be sufficiently large that an NLSP slepton would decay outside the detector. In 
this case, heavy charged particle pair production without missing energy could 
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be a manifestation of gauge-mediated low-energy SUSY-breaking models. 



Higgsino LSP 



The Higgsino LSP mode]ll23 involves a region of MSSM parameter space in 
which the X2 photino-like and the Xi is Higgsino-like, so the radiative deca^ 
X2 ~^ 7X1 dominates over other X2 decay modes (see, for example, Eq. (||)).lliZl 
The event can be again interpreted as (i) ee* production, but with e ex^, or 
(m) XiXi production, with Xi ~^ e~i?X2; ^^id the subsequent radiative decay 
of the X2 yielding the observed signature. 

In these models, photons only arise from the decay of X2- Other signatures 
involving two photons might come from the process i^v* — > 1^1^X2X2^ but there is 
no guarantee that the D is light enough to produce a substantial signal. Because 
the X2 is photino-like, direct X2X2 production is not large. In this model, the 
dominant neutralino and chargino production processes are XiXitXiXs: ^^'^ 
fxs- None of these involve the direct production of X2- Typically, the 



xtxi,x 



decay X3 Z*Xi occurs, yielding no photon. One of the next largest processes 
is xtx2j which would produce a trilepton signature in SUGRA models, but 
can produce i^lfix '^^ jjltx signatures in the Higgsino LSP model. 

If the stop is light, this discussion changes, because the Xi can decay 
Xi bii, followed by <i cxi- The signature is then a rather distinct 
"fhc$rp. However, such a light stop would appear in top decays, depleting the 
observed SM decays to an unacceptable level. This is only true, though, if 
there are no other sources of top quark production from SUSY, which there 
obviously can be. Surprisingly, such ptaodels can be constructed that are in 



accord with the present SUSY limitsl!-i3 If the gluinos are heavy enough so 



that g ^ tt\, or — > tti, and gluino production is further fed by squark decay 
Q — > qg, then one can compensate for the lost top quarks in SUSY decay 
modes. This leads to more sources of jbc$!rp events than just Xi^X2 events, as 
well as other signatures. 



Inclusive Two Photons and fJj. Signatures 

The generic "f"f$!x +X signature has no significant background from real pho- 
tons. The main backgrounds are caused by jets and electrons faking photons. 
The SM production of W{~^ ev)^ plus jets can fake some of the signatures if 
the electron is misidentified as a photon. These events have a Ij^rp spectrum 
typical of W events, peaked at about Afv^/Z ~ 40 GeV, with a long tail to 
high The dominant instrumental background, however, is from di-jet and 
7-l-jet production, where the large production cross section overcomes the small 
probability (~ 10~^ — 10"'^) that a jet fakes a photon. 
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Figure ^ s lasKS—the distributions from D0 (left) and CDF (right) 



diphoton events Eijilrj after imposing the selection criteria given in Table |16|. 
For the D0 analysis, the shape of the Iprp spectra agrees well with backgrounds 
containing two electromagnetic-like clusters, where at least one of the two 
clusters fails the photon selection criteria. Two events satisfy all selection 
criteria, with a predicted background, dominated by jets faking photons, of 
2.3 ±0.9 events. For the CDF analysis, the shape of the distribution is 
in good agreement with the resolution of the Z — > e^e~ control sample. The 
event on the tail in Ipj. is the "6677^2^" event. If the source of this event is 
an anomalously large WW'yy production cross section that yields one event in 
Mjjl^rp, CDF would expect dozens of events with two photons and four jets. 
However, the jet multiplicity spectrum in diphoton events is well-modeled by 
an exponential, and there are no diphoton events with 3 or 4 jets. As mentioned 
before, events with diphotons, jets and $!rp can be signatures of gauge-mediated 
low-energy supersymmetry-breaking models. 



Table 16: Selection criteria for searches 



Quantity 


D0 


CDF 




> 20, 12 GeV 


> 25, 25 GeV 




< 1.2 or between 1.5 and 2.0 


< 1.1 




> 25 GeV 


> 35 GeV 


A0 between 

If] J, and nearest jet 


N.A. 


> 10° 



D0 presents limitsE£2l in the framework of the Gravitino LSP scenario by 
considering neutralino and chargino pair production. Assuming M2 — 2Mi 
and large values of 'ti—, the signatures are a function of only M2, fJ-, and 
tan /3. Event rates are predicted using PYTHIA.Ei Figure |l^ shows the limit 
on the cross section for XiXi ^tnd X1X2 production as a function of the xt 
mass when is large and thus the xt mass is approximately twice the Xi 
mass. The figure also shows, more generally, the excluded region in the M2-/J, 
plane < gives larger xf 1X2 ~ Xi mass splittings, small |^| means xf jXij 
and X2 9.re more Higgsino-like) , along with a prediction for the region that 
might explain the CDF eejjfrp event as chargino pair production. The latter 
explanation requires 100 GeV < M~± < 150 GeV with M 



Xi 



113 



< 0.6M~± to 

Xi 



produce one event with a reasonable probability.!: 

As can be seen from Fig. |l^, the cross section limit is typically 0.24 pb for 
either XiXi or xf X2 production. By combining all chargino and neutralino 
pair production processes, a xt with mass below 150 GeV is excluded. Hence, 
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Figure 17: (Left) The spectra in the D0-search for events with 2 photons, one with 
Et> 25 GeV, the other with Et> 12 GeVE£3 The points are the data, the soUd line is 
the estimated background from di-jet events and direct photon events. The dotted Unes 
are for gaugino production within gauge— mediated models using the parameters listed and 
Ml ~ 2M2. (Right) The CDF ^j, spectrum for events with two central photons with 
> 25 GeV. Events which have any jet with > 10 GeV pointing within 10 degrees 
in azimuth of the l^j, are removed. The solid histogram shows the resolution from the 
Z — > e+e~ control saii*pl«. The dashed line shows the expected distribution from all SUSY 



production in a model tij with 



225 GeV, 
GeV. 



: 300 GeV, tan/3 = 1.5, and M 



Q 



300 



to keep the chargino interpretation of the eejj^J'p event, it is necessary to 



expand on the analysis of Ref. |113]. The chargino mass hniit is much higher 
than in SUGRA models, because of the 100% branching fraction for the de- 
cay Xi ~^ iG and the high detectability of the photon and ]}]r^. The result 
eliminates the possibility of observing signatures of this particular model at 
LEP200. The Iflrp cut needed to control QCD backgrounds makes the analysis 
sensitive to the mass splittings between Xi and Xi or X2- However, the sim- 
plest models predict unification mass relations between the gauginos, which 
thus gives acceptable mass splittings. 

D0 also has a limit on the cross section for ee* — > e~e~^X2X2i ~^ 
Z/1/X2X2J and X2X2 ~* 77X1X1 using the same analysis as for the Gravitino LSP 
search. Such signatures might also be expected in Higgsino LSP models. The 
limit on the cross section for such processes is about 0.35 pb for M~o — M~o > 
30 GeV, which is close to the maximum cross section predicted in these models. 
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Figure 18: (Loft) The D0 cross section limit on Xj^Xi ^iid X1X2 production, assuming 

M~± ?» 2Af~o and BR(x5' ^ iG) = 100%. The top dotted (dashed) curve is the cross 

^1 ^1 ^ ^ ^ ^ 

section from PYTHIA for X1X2 (XiX\) production. The bottom dotted (dashed) curve is the 
cross section Umit from the D0 collaboration|]£3 on XiX% {X\Xi) production. The vertical, 
hatched line marks the 95% C.L. lower limit on the lightest chargino mass from considering 
all chargino and neutralino pair production processes and all values of /i. (Right) The 
limits on the parameters M2 and in gauge— mediated models bassd-pn PYTHIA for tan /3=2 

and M~=800 GeV.t23 The hatched area is the region proposedtij to explain the CDF 
Q 

eeff^j, event. The solid line shows the D0 bounds. The long— dashed line shows a contour 
with M~+ = 150 GeV and the dash— dotted line shows a contour with M~„ = 75 GeV. The 

dotted lines show an interpretation of preliminary LEP results at an energy of 161 GeV. 



Single Photon, Heavy Flavor, and fj 



the signature jbcJ^j,^ as predicted in Higgsino LSP mod- 
The data sample of 85 pb~^ contains events with an 



CDF has searched fo 
els with a light stop 
isolated photon with > 25 GeV and a jet with an SVX 6-tag. The spec- 
trum of these JiiZEnts can be seen in Fig. 19, After requiring ^j,>20 GeV, 98 
events remain 



11£ 



The estimated background to the 98 events is 77 ± 23 ± 20 events. The 
shape is consistent with background. About 60% of the background is due to 
jets faking photons, 13% to real photons and fake &-tags, and the remainder 
to SM jbb and jcc production; all of these sources require fake ]^rp. When 
the cut is increased to 40 GeV, 2 events remain. More than 6.43 events of 
anomalous production in this topology is excluded. 

The efficiency used in the limits is derived from a "baseline" model with 



M- 



= 40 GeV, M~o 

Ai A2 



70 GeV, TTif^ = 60 GeV, m~ = 250 GeV, and Mg 
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Table 17: Summary of the 85 pb ^ data sample for the CDF 'yb'$j< search. Limits are set 



using all cuts which results in two events. 



Quantity 


Cut 


Cumulative Number of Events 


£;J, ID cuts 


> 25 GeV 


511335 


SVX 5-tag 


> 1 


1487 


E!^,, \r]\ < 2.0 


> 30 GcV 


1175 




> 20 GeV 


98 


A(^(7-^-r) < 2.93 


> 40 GeV 


2 



225 GeV.0 The distribution of the number of jets in the data is shown in 
Fig. |l^ compared to that expected from backgrounds and the SUSY model 
(scaled xlO). There are more jets expected in the SUSY model than the data 
indicates because of the hard kinematics of squark and giuino decays. The 
baseline model predicts 6.65 events, so this model is excluded (at the 95% 
C.L.). This result does not rule out the Higgsino LSP model in general, only 
one version with a fairly light mass spectrum. A more general limit can be 
set by holding the lighter sparticle masses constant and varying the squark 
and giuino masses. In this case squarks and gluinos less than 200 GeV and 
225 GeV, respectively, are excluded. 

4. 10 Other Anomalies 

There are other anomalies in the current data beyond the "6677^^^" event. 
These are, so far, either single, rare events or discrepancies on the tails of dis- 
tributions where statistics are low and backgrounds difficult to calculate. In 
addition, there is the problem of calculating probabilities for "anomalies" a 
posteriori. The expected number of events in any one channel from SUSY is 
usually small with the present integrated luminosity. New physics will most 
likely show up as a few events on the tails of SM distributions. Since there 
are many potential SUSY signatures, one can only follow a strategy of sys- 
tematically analyzing all high-mass channels and looking for discrepancies on 
the tails of distributions. The single events such as the "6677^^-" event have 
been useful as "guideposts" indicating promising new channels, such as the 
^hj$j. channel described above. It is still possible that a sensible picture of 
these events will emerge from the Run I data when a complete survey of all 
channels is completed using both detectors. At the very least, this is an im- 
portant exercise for preparing the Run II analyses. 

"This analysis predates LEP results which exclude this example. 
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Figure 19: CDF results for the $j, (left) and the jot multiplicity with Et > 15 GeV (right) 
in events with a photon and a SVX 6-tag. The search is for the signature b-yfrp in a scenario 
where X2 ~' 7Xi ^^^d the stop is hght. The jet multiplicity histogram is made by requiring 
-^y>20 GeV. The SUSY model is normalized to the area of the data histogram - this is 
scaling by a factor of 100 for the histogram fttwi a factor of 10 for the Jijot histogram. 
The SUSY model has a Higgsino LSP 



generated with PYTHIA 6.1. 



Top Dilepton Events 

As discussed earlier, the signature of dileptons+2 iets+$!j^ is a promising SUSY 
search channel (see Sec. ^2|). However, such events would also be a background 
to the SM top quark search using dileptons. The consistency of this dilepton 
sample with tin at expected from tt production has been the subject of intense 
investigationlHil There are a number of peculiarities, none by themselves sta- 
tistically significant at a level required to claim new physics. However, there 
are several eventsjthat have low probabilities of being from top decay or any 
other SM processll£3 Such events should be taken seriously as potential SUSY 
candidates. 

The most interesting of the anomalous CDF events 
of Run 67581, which has three, clean, isolated, high-pT^ leptons, large $Jj,, and 
a high-i?T jet. In addition, the most energetic of the leptons is a positron with 
Et — 200 GeV, significantly larger than is typical for top events (0.06 ± 0.02 
events are expected). The corrected is over 100 GeV, also large for top 
decay ( 0.6 ±0.1 events are expected). The event contains a jet with Et — 
100 GeV; the total transverse energy plus l^j, is about 450 GeV. The other 
two leptons are an electron with Et — 27 GeV and a muon (/x~) with pT = 27 



53 is Event 129896 
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GeV/c. The invariant mass of the e^e^ pair is 130 GeV, well away from Mz; 
the pair has very high pT- In the SM top quark analysis, the event is classified 
as a dilepton+2 jet event, because the lower Et electron fails the fiducial cut 
by 4 mm and is thus defined to be a jet; however, the electron passes all other 
standard electron criteria and is a "golden" electron in all other ways.^ The 
kinematics of the event are unusual: the invariant mass Me^j is on the order of 
nit , while the other hemisphere contains only the high-ET positron. The three 
isolated leptons and the kinematics make the event unlikely to come from SM 
top production and decay. The event is aJiigh-mass trilcpton +frp event, and 
is consequently a good SUSY candidateE 



data: ★ e/i O ee □ 
tf MC 

• Z->ttMC 



CDF (1 10 pb-') 






*lit, SLT-SLT 








eT,SECVTX 

Cut 








candidates after all cuts 




data otter lepton ID cuts only- 







[GeV] 0, IGoVl 

Figure 20: (Left) Scatterplot of the angle A<j){f!rp,ij) between the corrected and the 
closest lepton or jet versus corrected for the ee, fifi, and efi candidate events, compared 
with the expected distributions for tt and background. Background and top contributions are 
not normalized to the expected number of events. (Right) The distribution of $rp significance 
versus ^j, for events with a primary lepton and a tau candidate (the slope of the data and 
background is different, because the background is dominated by QCD) in the CDF data 
compared with the tt Monte Carlo. Three of the four final candidate events (stars) have 

fe-tagged jets. 



Other discrepancies in the top dilepton sample involve the kinematics. 
Some of the anomalous behavior in the kinematics can be seen in Fig. Heft), 
which shows $!rp versus A<j) between the l^j, and the nearest jet or leptonHllH^ 
Also shown is the distribution expected from Monte Carlo tt events, but corre- 

°In the top quark analysis the fiducial volume was conservatively chosen to be the same 
as for the precision electroweak measurement of the ratio oi W to Z cross sections. The 4 
mm miss does not affect the electron identification. 
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spending to 100 times the luminosity. There are several events out in regions 



less populated by top quark events (one is the tdlepton event). Figure 20 
(right) shows the distribution in Jpj, significance E3 vs Iprp for the CDF tau- 
lepton top sample. None of these latter discrepancies is at a statistical level 
that is significant; these will be channels of great interest in Run II. 



5 Conclusions 

As can be seen from Table 0, there has been a large effort in SUSY searches at 
the Tevatron. However, given the wide range of possible experimental signa- 
tures in the minimal Supersymmetric extension of the Standard Model, there 
is still work in progress and much to be done. Many Run I analyses are under 
way. 

Our quantitative conclusions on Run I are reflected in the Figures and 
Tables of this review; here we will add a few more general qualitative observa- 
tions: 

1. A systematic exploration of signatures and channels is just starting. In 
addition, the detectors have not yet been exploited fully; for example, 
better c~tagging and dijet resolution to investigate flnal states with re- 
constructed W and Z bosons may be possible. These tools will allow the 
study of new channels. 

2. There are some events involving leptons and/or photons that are provoca- 
tive, and can be "guideposts" for Run II and further Run I analyses. 

3. There is a substantial need for theorists and experimentalists to work 
together to understand better how to derive and present limits from 
the Tevatron. We should move on parallel paths toward more "model- 
independent" predictions and limits {e.g. presenting plots of cross section 
versus experimentally measured quantities like thresholds), and confront 
specific models in ways that allow the two experiments to compare their 
results. 



4. The analyses so far are luminosity limited: the reach of the searches is 
just entering the interesting regions. 
In Run II, two upgraded detectors at the Tevatron will collect more data 
at a higher energy of 2 TeV. The nominal integrated luminosity is 2 fb~^, with 
a possible extension to 10 or even 30 fb^^. The production cross sections for 
heavy sparticles will increase significantly with the higher energy. Chargino 
and neutralino searches, as well as squark and gluino searches, will cover a wide 
range of SUSY parameter space in Run II. Most importantly, by extending 
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Run II up to an integrated luminosity of about 20 fb^^ and combining search 
channels, the Tevatron can perform a crucial test of the MSSM Higgs boson 
sector. 

The experience gairied-frnm Run I analyses will greatly increase the quality 
of the Run II searches. New triggering capabilities will open up previ- 



ously inaccessible channels, particularly those involving r's and heavy flavor. 
Increased 6-tagging efficiency and If^rp resolution will enhance many analyses. 
A factor of 20 or more data combined with improved detector capabilities 
makes the next Run at the Tevatron an exciting prospect. 
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A Appendix A 

A.l Typical Decay Modes of Supersymmetric Particles. 



Table 18: Typical final states from sparticlc decay, assuming x">X^>^>'^ < ©{t^ *>^)i9- 
HLSP denotes models with a Higgsino LSP and GLSP denotes models with a Gravitino 
LSP. Event signatures from sparticle pair production can be constructed by combining two 

decays. 



Particle 


Intermediate State 


Final State 


Comment 




ti* 
^bb* 


$T 

jj $T 

bWc $j. 

bWM $T 
^ bWbjj $r 
bb $rp 


HLSP, GLSP 
mi<m\,-\- M~± 


xt 


bi* 
-^tb* 


331 
be $T 

bU $T 
^ hbjj $T 
bWb 


HLSP, GLSP 
HLSP, GLSP 
nif <mb + M~± 


e 


^^X°2 

^xt 


id 

^ ^33 $T 

^ 33 $T 
be $rj. 


HLSP, GLSP 

rui <mb + M~± 
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A. 2 Typical Decay Modes of Supersymmetric Particles (continued,!). 



Tal>ic 19: Typical iiiial sla,lcs iix.)iu b])arUck' decays (c-oiiliiiuecl.l). 



Particle 


Intermediate State 


Final State 


Comment 






(7 $T 








^ 3 J $T 

el $j. 

£bc $T 


HLSP, GLSP 

mi <mb + M~± 


t 






mi <m.b + M~± 




^bxt 


bl $T 

^ bj] $T 
bW $T 


Ml >mt + M-o 




^bW 


^bW $j. 




b 


^bxl 
^bx^ 

'^txf 


~>b $j, 

^b $T 
bU $T 

bjj $T 
bWi 

bWjj $T 
bWbc 


mi <m,b + M~± 






cW $T 
b£W $T 

^ bjjW $T 


mj < rofe + M~± 

Xi 
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A. 3 Typical Decay Modes of Supersymmetric Particles (continued, 2). 



Tal>i(.! 20: Tyi)i(.*al liiial sla,lcs i'roiu bi)arUck' decays (c-oiiliiiuecl.2). 



Particle 


Intermediate State 


Final State 


Comment 


Q 










^3X1 


^3$T 
















333 $T 






^3xt 


^ 3^ $T 








333 tr 








3 be $T 


mi<mh + M~± 




ifi 


^ J33 $T 








— 3' I t'v 


IILSP. GLSP 


9 


^jQ 

tt* 

bb* 

-> 33Xi 
33X1 

33xt 

^ ttXi 
^tbxf 


^ 33 $T 
bWc $T 

bb $j, 

^ 33 $T 

^ 33 $T 

-^33U$T 
3333 $T 

^ 33^ $T 
3333 $T 
jjbc $T 
bWbW $T 
bWbe $j. 

^ hWbjj $T 


TO£ <mb + M~± 
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B Appendix B 

B.l Examples of Multijet Signatures for SUSY. 



Table 21: Examples of R— Parity Conserving SUSY signatures at the Tevatron: Jets+-^jn. 

Not all signatures arc listed - we have (somewhat arbitrarily) restricted the list. We assume 
that the LSP is the Xi- Note that Q decay's give one or 3 jets. (/ decay's give 2.1. or 6 jets, 
aiid the x'l^ ■ 'Hi^l dcca.\'s gi\x; a,ii even iiuiiiljt;r of jels. 



R-Parity Conserving Signatures: Jets 


Signature 


Production 


Decay 


3$T 


Qx\ 


Q^9X? 


jj$T 


QQ* 


Q ^ 9X?; Q* ^ 9X? 




xfx? 


xt -> mxi 




ti* 




jjjfr 


Qx? 


Q 9xf ,xf QQXi 




Q~9 


Q ^ qxi;g qqxi 




QQ* 


Q ^ 9X?; Q* ^ qxt ^ q{qqx°i) 




99 


g_-^ qqxi 




Qg 


Q qxt q{qqxi)\ 9 qqxi 




QQ* 


Q qxf q{qqxi) 




99 


9^-^ qqxi; g ^ qQ ^ 9(gxf ^ q{qqxi)) 


> &3$T 


Qg 


Q qxt 9(99x1); 9^ qQ-> q{qxt ^ 9(99x1)) 




99 


g^ti^ iWb){xtb) - {jjb){jjbx1y, ~9 - 99X? 
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B.2 Examples of SUSY Signatures that Include b-quarks. 



Tabic 22: Examples of R— Parity Conserving SUSY signatures at the Tevatron: b-tags + jets 
We have shown only a few modes. The signatures change depending on the relative 
masses of the 6, i, xt ^-nd X? ■ 



R-Parity Conserving Signatures: b quarks 


Signature 


Production 


Decay 


b^T 

bjji^T 
bjjj . . 


It 

bg 

bb 
bi 
bb 


b^bxl 

b^bx'i;i^cx'i 
b->bx";g^ qqXi 
i^bxt 

b^bxl; b~^iW,i^cxi 
b^bg,g ^ qqXilJ ^ cx? 
b^bgg^ qqxl; b iWi cXi 


bbj ...fj- 
bbbj . . 
bbbbj . . .$T 


bb 
tt* 
fb 

gi 
gg 


b^bxl 

t^bxf,xt^qq'x°i 

g bbx\\ bxi 

g bbxi; i^ bxt,xt qq'xi 
~g - bbx°. 


bbljj$j. 
bbU$T 


tt 

a 


i ~^ ' ^ ej/x^/, xf qq'Xi 
i^bxi,xt ei^xi 
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B.3 Examples of SUSY Signatures that Include Leptons. 



Table 23: Examples of R— Parity Conserving SUSY signatures at the Tevatron: Leptons. We 
have shown only a few modes. Note that (for example) e decays can give 0, 1, or 3 (charged) 
leptons, and decays can give 1 or 2 leptons. The signatures will change depending on the 
relative masses of the sneutrinos and sleptons in the different generations. We have not shown 
explicitly the differences in the "left" and "right" slepton decays. The decays that single out 
the T (for example, from the H'^) are omitted here. Decay modes involving neutralinos and 
charginos can be created by feed-down from squark and gluino decays. Gluino decays can 
lead to leptons with uncorrclatcd charges. 



R-Parity Conserving Signatures - "Generic" Leptons 



Signature 



1^ 



m...3...$j 



Production 



ev 

113 

Xi xt 
xlxl 

ee 

xtxl 
xlxl 

ee 
ev 



xfx" 

ee 

xfx2 

~99 

it* 

xfxs 
99 



Decay 



Xi 

x^ 



xt 

x^ 



ei/xi; X2 

evXi _ 

eex?; X2 

ex? ^ 

evxi] xl 

eexi; X2 
-0 ~o 



6X2^X2 eexi 



v^Xi 



eexi 

eexi 




w 

9^' 
e — s 

xt 

9^' 
i 

xt 



eX2>X2 



eexi; i> 



ext 



^75 — -~o 

-» euxA_; X2 QQXi 

' Q<^xt'^ 9 ^mxi 

' i ^^'^ti^^^xi) 

99'Xij^ X2 eexi 
• «;9'xf , xt ei/Xi 
bxt^xt ei/Xi ^ 
^ ej^xi; X3 mxli-^ eexi) 



m'xt'-: 9 



qqx2 
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B.4 Examples of SUSY Signatures that Include Photons. 



Table 24: Examples of R-Parity Conserving SUSY signatures at the Tevatron: j's+flrp. 



In the light Gravitino scenario (GLSP) x? ^ "iG always occurs. If the decay has a long 
lifetime, one of the two Xi may decay outside the detector (LLC). In the Higgsino LSP 
scenario (HLSP), X2 ~* 7Xi often occurs. 



R-Parity Conserving Signatures: Photons 


Signature 


Production 


Decay 


Comment 


1333 ■ ■ -^T 


Qxl 
xtxl 

~± ~o 
Xi X2 

99 

Q~9 


X?^7G,x?^7G 

Q^gx;.x?^7G 

xt ^Wx\,x\^lG 

xt ^m'x\,xl 1x1 

g-^ qQ q{qxt ^ qilQXi)) 

Q gx?; 9 qqxl,xl ixi 


LLG 

LLC 
LLG 
HLSP 
LLG 

HLSP 


lb$T 


~± ~o 
Xi X2 


Xf ^ th,t ^ cxl; X2^7Xi 


HLSP 


1^$T 


~± ~o 
Xi Xi 
~± ~o 

eu 
99 


xt evxi,Xi iG 

xt ei^Xi; X2 1x1 

e ex2, X2 7X1; vxi 

9 QQ'Xi{-^ eiyxi); 9 ^ 99X2' X2 ^ 7Xi 


LLG 

HLSP 
HLSP 
HLSP 




X2X2 
99 


X2 ^ eexi; X2 ^ 7Xi 

9 99X2 > X2 eexi, X2 1X1 


HLSP 
HLSP 


1133 ■ ■ -tr 
llijj ■ ■ -Vt 
-yyUjj . . .$T 


X?X2 

~it~0 

Xi xi 

ee 

xtxi 


xt qq'xi; X2 qqxi, x? ^ iG 

xt euxi; X2 ^ QQXi, Xi ^ iG 
e exi,Xi iG 
xt ez/X2,X2 ^ 7X1 


GLSP 
GLSP 
GLSP 
HLSP 
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